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ABSTRACT 


This thesis is devoted to the study of mathematical models 
which have been proposed to describe the activities of the nerve cell. 


The models considered can be classified into two basic categories. 


The models in the first category are those which describe the 
electrical behavior of the nerve membrane. Several mathematical model 
which describe the electrical behavior of the Loligo giant nerve axon 
membrane are considered. It will be shown that the models satisfy a 
generalized formulation. These models for nerve excitation, impulse 
initiation, are extended to incorporate the phenomena of impulse 
propagation in the nerve. One of the excitation models, that proposed 
by Hodgkin and Huxley (1952d),is modified to describe excitation and 
propagation in myelinated nerves. ‘Pinally two models are considered which 
account for the affect of the geometry of the dendrites on synaptic 


input. 


The second category of mathematical models are those which 
describe the repetitive impulse activity of the nerve. Three such models 
are considered. As in the case of the excitation models these will be 
shown to be special cases of a generalized model. By analysing this 
general model it is possible to determine an interaction equation for 
the activity of nerves in a nerve network. This mathematical nerve 
network model forms the basis for three nerve network models which 
are considered. For one of these, the linear behavior of the solutions 


is examined for small nerve networks. 
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INTRODUCTION 


The purpose of this thesis is to survey the mathematical models 
which have been proposed to describe the activities of a single nerve 


cell. 


The fundamental characteristic of the nerve cell is the nerve 
impulse. Functionally and morphologically the nerve itself can be 
divided into three parts; the axon, the soma, and the dendrites. A 
nerve impulse is generated in the axon portion of the nerve. The axon 
of one nerve impinges on the dendrites and soma of other nerves at 
points called synapses. Once an impulse has been generated in the axon, 
it will propagate over the entire axon and its effect will be transmitted 
to the soma-dendritic region of the other nerves through these synapses. 
The dendrites and soma of the nerve therefore act as a receptive field 
for stimulation of the nerve. In vitro, if sufficient stimulation is 
received in the soma-dendritic portion of the nerve, an impulse is 


generated in the axon portion. 


The nerve impulse is essentially an electrical phenomena which 
can be recorded using microelectrodes. Furthermore, stimulation of 
the nerve by electrical current of a sufficient strength will initiate 


an impulse. 


The physiological properties of the nerve originate in the 
membrane which forms its surface. Using the techniques of intracellular 


recording, extracellular recording and voltage-clamp it is possible to 
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determine the electrical properties of the nerve membrane. Chapter I 
of the thesis is devoted to a survey of mathematical models which have 
been formulated to generate this electrical behavior, in particular 


the initiation and propagation of a nerve impulse. 


In certain cases the only information which can be obtained 
regarding the nerve is the pattern of impulses produced. Such 
impulse patterns show a great deal of variability in the interval 
between impulses. In Chapter II, mathematical models are considered which 
have been constructed to generate trains of impulses that resemble the 


experimentally observed nerve activity. 


In studying the properties of the membrane and the impulse 
patterns, a wealth of experimental data has been assembled. The basic 
objectives of the mathematical formulations are to provide adequate 
theoretical concepts for the processes involved in the production of 
such behavior as well as methods for the analysis and comparison of 
this data. Using these derived formulations it is possible to predict 
the behavior of the proposed model, in this way providing a means of 


testing the model and establishing guidelines for further experimentation. 
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CHAPTER I 


Introduction 


In Chapter I the mathematical models which have been formulated 
to describe the electrical behavior of the nerve membrane are discussed. 
As stated in the Introduction, the initiation of the nerve impulse is 
associated essentially with the axon portion of the nerve. Since the 
impulse is an electrical phenomena, it is therefore ia eens to 
consider the behavior of the axon membrane to electrical stimulation. 

The axon considered will be Loligo giant nerve axon which has been 


extensively studied. 


In the absence of stimulation, there is a constant potential 
difference maintained across the membrane. If a stimulus in the form 
of a brief outward current pulse is applied to the nerve axon through a 
stimulating electrode touching the membrane, the recorded membrane 
potential changes from that of the resting potential. The magnitude 
of this change is dependent upon the amplitude of the applied stimulus. 
For a very weak stimulation, there is a temporary change of potential 
which decays away not only with time but also with distance from the 
electrode. With increasing stimulus strength the recorded potential 
change increases faster than the stimulus amplitude. The resulting 
potentials still decay both temporally and spacially. Finally a critical 
level of stimulus amplitude, called the threshold, is reached. If a 


stimulation, with this amplitude is administered to the axon, the recorded 
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potential curve increases abruptly to form an unmistakable triangular 
wave, the action potential, which does not decay with time and distance. 
For stimulus amplitudesabove this level, which are not biologically 
damaging to the tissue, an action potential results. The action potential, a 
nerve impulse, is initiated near the electrode but immediately splits 
into two waves which travel at a constant waveform and velocity in 
opposite directions away from the electrode over the nerve. For 

natural stimulation through the soma-dendritic region, the action poten- 
tial is a single wave which propagates over the entire nerve. Thus the 
action potential once initiated can be recorded at any point on the 
nerve as it passes. It is however, acee ie to determine from such 

a recording where it has originated on the axon or the strength of the 
applied stimulus. Nerve excitation, the initiation of a nerve impulse, 
therefore is an all-or-none phenomena for which only the presence of 


absence without intermediates, can be deduced. 


If a second above threshold stimulation is administered during 
a short time interval immediately following the action potential, it 
is possible to initiate another impulse. This time interval is called 
the absolute refractory period. Following this is a period of time during 
which an action potential can be initiated by the second stimulation but 
the threshold value in this case is above that for the resting nerve. 
During this relative refractory period the threshold level gradually 


returns to its resting state value. 


If the first stimulus is below threshold for excitation, a 


nerve impulse is not produced, however, the membrane is affected. Ina 


an Sis LN, Dm fe hisae’ 


nelogesisd sldstateiauy as weed 02 sci eeesezont svr03 J bata 
sodeseio bee sents Aste saa den, ses dokide , setIno70G' —n fa- 
¢iisskgolerd ton ars doide (favel etrs avoda eeiad Eben | 


Tey, te ve 
fattnsidg Horsos, Sit .2 Lips winivmaetog nolsos fie saueeks |ad2 oF ie 


eaiiqa ¢tisisibeans tud sboxsisibe act +son S4tnidihdyeb nue 
nk vos soley. bie ctoboviee auadentts a te Cayeal “rho bithe aavew Os a 

yOu -~svxse snd wsve ehorasess ait mots dock (anekaoeniel z 
~wsitoer motions sit. ,Antaes atte bdebuanbe “sd? deinen hokse tints 
Sis aaiT .oyyst suitng Sarto Soyo bay aaeueys Deki SVEW signta ee 


J _ 
ime ) 


BAF 630199 Uns sae babrbgo% 3d aid hana site wore feb a 
ope most esiuseieb ‘oF iraisdaaa ,yavsvod, af FE t -eaeaEg: ail 
afi) to\Adguesd8 $2) 26. coxs sd) Bo bas ack atac: wad 3£' disoite 

 seivqut avien & 20 moisetatat sad eat te: pvaet, -eudeiniy2 |f : 
to, bonseeye SH2 yino singe 1203 Brtaadnadt Sabin 25- Lo, os 5 ie 


bsonbab 93° ne /uoasthearzaaat Suodan 


an 
i is Pash si 


etitws bsysselatmbs: at Rot Euan Se bisdes sii wis: bugsee ‘S ay 
rf Sebieatog obses oid agutivol tor yisgabboumet, lavtonat wm 
beLbes et Levaaaak emt ati?) sesiuqah satiroae alaldank oF, 
cate id 

arith sata: io ities BL) abet gatwottos pbokreg igrofoaviss — 
ind rods tums hnooee. ett vw peyana int ad Bo ted ge30q) holden: ain 


aimee 108: oe avods ab sac adaduenil 


short period following this subthreshold stimulation, the threshold 
value for a second stimulation will be lower than the resting value. 
After this phase, is a period of time during which the threshold for 


a second stimulus is greater than the resting value. 


A stimulus applied to the nerve therefore has a dual effect. 
The first is excitatory, which results in a propagating nerve impulse 
if the stimulus is strong enough. Following this, the effect of the 


stimulus is to reduce the excitability of the nerve. 


The mathematical models of nerve excitation are formulated 
in an effort to determine the variables and interconnections which 
result such electrical behavior. The major objective is to determine 
the construction and operation of the membrane. As shall be seen it 
has been possible in certain cases to associate with the mathematical 
model, a physical or chemical model of the membrane which satisfies 
the proposed mathematical description. Although the major objective 
has not been attained, the models have provided guidelines for 
physiological research, and have assisted in translating physiological 


concepts into quantitative and logical terms. 
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§I.1 Early Mathematical Models for Nerve Excitation 


The first mathematical models were formulated in an attempt to 
explain the observed behavior of the nerve to an applied electrical 
stimulus. Experimentation made it possible to determine a number of 
quantitative relationships between the variables of the stimulus and the 
initiation of an action potential in the nerve. It was these relation- 


ships which provided the basis on which to test the models. 


An early model which was extensively examined for its 
predictions regarding the strength-duration relationship in the nerve 
was proposed by Blair (1932). Blair postulated that the initiation of 
an action potential by an applied electrical current was governed by an 
"excitation toedenes in the nerve which was under the influence of the 
stimulus. The state of this process was assumed to approach a threshold 
level for excitation at a rate proportional to the applied voltage; and 
if the stimulus was removed before the threshold was reached, to decrease 
towards a resting level at a rate proportional to the state of the pro- 
cess. Thus the dynamic behavior of the state K of the "excitation 


process" satisfied 


(i142) & = bv - kK 


where V is the applied voltage, t is time, K=0 is the resting 
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level, and both b and k are constants. An action potential was assumed 
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The formulation given by (1.1.1) and (1.1.2) had originally 
been considered as a model for nerve excitation by Lapique (1907). The 
difficiency of the model was in its inaccurate formulation of the 


strength-duration relationship for a step current. 


The solution of (1.1.1) for stimulation by a step current 
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if the "excitation process" is at the resting level at t=O. By 
setting K =L in (1.1.3), the strength-duration relationship for a 


step current is given by 


= kL 
V 


where R is the rheobase voltage. The major difference between 


O 
the computed graph of log cents) considered as a function of t_ for 
Oo 


various values of k and the experimentally obtained graphs is that 


the curves in the latter case do not all pass through the origin. These 
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results suggested to Blair that the strength-duration relationship for 


a step current is of the form 


V 
log 7 =kt + C 


where C is a constant. 


To theoretically justify the strength-duration relationship, 


Blair proposed the threshold L must be a function of V_ given by 
Gia ls4) L=a + BV 


where Q and 8 are constants. 


This revised formulation consisting of (I.1.1) and (1.1.4) 
permitted good approximations of the experimentally obtained strength- 
duration curves for stimulation by a step current, a linearly rising 
current with a sufficient rising rate, and condenser discharge (Blair: 
1932-35). Unfortunately the theory was unable to account for a number 
of electrophysiological facts; in particular, the equations did not 
satisfy the constant quantity relationship for rectangular current of 
short duration (Hill: 1936b), and linearly rising current would always 
result in excitation irrespective of the gradient of the rise (Rashevsky: 
1938). Furthermore, the assumption that the change of threshold depends 
only on the instantaneous voltage leads to conclusions which are con- 
trary to known results (Rushton: 1934). These difficulties, Rushton 
suggested, result from the fact ene theory does not adequately account 


for the accommodation property of the nerve. 
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The difficulties were overcome by a model proposed by Hill 
(1936a). For this model, that property of the nerve in the neighborhood 
of the stimulating electrodes which was influenced by the electric 
current was called the "local potential". The ene behavior of the 
"local potential" with only slight modification was the same as that of 
the state of the “excitation process" in Blair's formulation. The 


"local potential" K was determined by 


Gal.) ee ae 


where I is the applied current, Ko is the resting level of the 
"local potential", and both b and k are constants. As in the model 
proposed by Blair, excitation was assumed to occur if the "local 


potential" reached the level L. 


The major difference between the two theories was in the for- 
mulation for the behavior of the threshold to an applied electrical sti- 
mulus. Blair assumed the threshold was determined by (1.1.4) while Hill 


proposed that 
(TL) a Rac ea Se 


where 8 and } are constant. The formulation given by (1.1.6) was 
adopted by Hill to account for accommodation in the nerve and the 
observation that the threshold gradually returned to a resting level 

Ly if the stimulus was removed. Hill simplified (1.1.6) by considering 
the particular case for which a nerve, described to be fully accommodated, 


responds as if no current were passing through it; that is, 
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K-L=K - Lo . Hence (1.1.6) reduces to 
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Hill's model given by (1.1.5) and (I.1.7) was shown to be 
sufficient to explain anode break excitation, and the necessity of a 
minimum current gradient for excitation by a linear rising current 
(Hill: 1936a). The theory gave satisfactory descriptions of strength- 
duration curves for a variety of stimulating current, and good approxi- 


mations of the results for numerous related experiments (Katz: 1939). 


The theory of Hill has been classified as a '"Two-Factor Theory" 
for nerve excitation. Hill's model is only one member of this class of 


models. 


Rashevsky (1933) also proposed a two-factor theory to explain 
the initiation of an action potential by an applied electrical stimulus. 


The mathematical formulation for the theory was 


AG =bI - k(K-K.) 
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where I is the stimulating current, and a, b, k and m are constants. 
In this theory K was chosen to represent the concentration of an 
excitatory ion in the neighborhood of the stimulus while L represented 
the concentration of an inhibitory ion near the stimulus. Ko and Lo 


represented the resting level concentration of the respective ions. 
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Excitation was assumed to occur if K>L. 


The formal difference between Hill's theory and Rashevsky's 
theory is that in the latter, the term which accounts for accommodation 
is a function of the applied current, while in the former the term is 


an implicit function of the applied stimulus. 


In 1934, Monnier proposed a two-factor theory for excitation 
which depended on a single process, the "state of excitation". During 
the application of a constant current I , the "state of excitation" 


was defined by 


CI. 229) . e = KI(e -e ) 
where K, Ty 5 and T, are constants. Hill (1936a) succeeded in 
showing that e« corresponded to the quantity (Ly ~ KD) - (L- K) of 


his theory, and that theoretical predictions regarding excitation for 


the two are equivalent. 


The final model to be considered in this section is the 
generalized model of the '"Two-Factor Theory" class proposed by Young 


(1937). The mathematical formulation was 
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where k k k k and b are constants, K represents the 
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excitatory process with resting level Ko » L represents the threshold 
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with resting level Ly » and I is the current strength. Young showed 
that, in terms of time-intensity relationships, Rashevsky's and Hill's 


theories were equivalent to this more general formulation. 


Although the two factor theories were a satisfactory repre- 
sentation of events leading to excitation by an electrical stimulation, 
they were unable, except in one case which will be considered in 81.5, 
to describe events following the initiation of the impulse. A further 
drawback of these theories is the fact that no physical basis existed 
for the variables which were chosen only to represent biological 
phenomena that had been experimentally observed. To overcome these 
difficulties more extensive knowledge of the nerve was required. This 
necessitated the advancement of experimental techniques and equipment. 
The most important development was the "voltage-clamp" technique of 
studying the nerve. This advancement lead to the single most important 
development in modelling nerve excitation, which is the subject of 
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81.2 The Hodgkin-Huxley Theory for Nerve Excitation 


The most significant development in describing nerve excitation 
and events associated with it, was obtained by Hodgkin and Huxley (1952 


a-d), initially working with Katz (Hodgkin, Huxley and Katz: 1952). 


Experimental results obtained by employing the "voltage-clamp" 
technique on the axon of the giant nerve fiber of "Loligo" indicated 
that a current flowing through the membrane of the axon could be 
separated into a capacitance and an ionic current. The ionic current 
consisted of three components: a potassium current 1. » a sodium 


current Ia , and a "leakage" current I) which was carried by chloride 


and other physically unidentified ions. 


Hodgkin and Huxley postulated the ionic current to be a result 
of the movement of the ions down their respective electrochemical 
potential gradients. Hence, the current resulting from each component 
could be defined in terms of the ease with which the component ions 
crossed the membrane, the permeability of the ions, and the difference 
between the membrane potential and the equilibrium potential for the 
ions. Since the instantaneous current-potential difference relationship 
for each component was found to be linear, the mathematical formulations 


describing the current densities of the components were 
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Ey and Eva are the equilibrium potentials for potassium and sodium, 
respectively, and x and na denote the conductance per unit area of 
membrane of the respective ions. Eo is the membrane potential at which 
the "leakage current" is zero and Bo denotes the conductance per unit 
area of membrane of the ions contributing to the "leakage current". The 


conductance term in each case is a measure of the permeability of the 


component ions. 


For practical application (1.2.1) was transformed into 
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and Ee is the resting membrane potential. 
Thus the total current density I per unit area of membrane is 


2 etal z 
CE.2.4) I= Cy ae + Buy, (V v0 + g(V - Vv) + 8, (V ~ Vo) 


where accounts for the capacitance current with C the measure- 


dv 
Onde. M 


ment of membrane capacitance per unit area, and t is time. 


V. 


The experimental results indicated C,, V Na? 


mM? *K? gust ek» 


are constant but and Byq are functions of both time and membrane 


eK 
potential. From the data Hodgkin and Huxley determined that for a given 
membrane potential the experimental curve of Be = g(t) was best 
described by assuming potassium conductance was proportional to the fourth 


power of a variable n which itself obeyed a first order differential 


equation. Hence 


= a 
By = By 
(ie2e5) 
dn = — -_ 
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where a is a constant which represents the maximum potassium con- 
ductance, and 0 Sen Sk a, and B are rate constants which were 
determined to be functions of the membrane potential for a fixed tem- 
perature and calcium concentration. In the case of sodium conductance 
it was proposed that the experimental curve of Sxa ™ Bua 6b) for a 
fixed membrane potential could be described best by assuming that Sua 


was a function of two variables m and h , both of which obeyed first 
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order differential equations. The proposed equation are 


Ke, 3 
eNa a Ena a 


| dm 
1.2.6 J aaa tS _ —_ 
( ) Was a (1 - m) - 6m 
| dh yy 
tl dt eit _) vee 
where Bna is a constant representing the maximum conductance, 


O<m<1 and O<h<1. As in the case of potassium conductance, 
for a fixed temperature and calcium concentration, the rate constants 
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me h and BF were found to be functions of the membrane potential. 


The functional relationships between the rate constants, in 
both (1.2.5) and (1.2.6), and the membrane potential were empirically 
determined to fit the experimental data and at the same time conform as 
closely as possible to the equation derived by Goldman (1943) for the 
movement of charged particles in a constant field. The resulting 


functions were 
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The reason for trying to conform to Goldman's equation which is of the 
form 
ae - d 
e - 1 
where C, a, b , and d are constants, will become clearer in the 


discussion of the proposed mechanism for membrane permeability changes. 


Thus the set of equations (1.2.4) - (1.2.12) makes it possible 
to determine the total membrane current I as a function of time and 
membrane potential. These equations are referred to in the literature 


as the "Hodgkin-Huxley equations". 


Hodgkin and Huxley, in addition to deriving an empirical 
formulation determining the membrane current, suggested a possible 
physical model to describe the potassium and sodium conductance changes 
in terms of the mathematical formulation. Potassium permeability was 
hypothesized to be controlled by potassium specific charged particles in 


the membrane, the movement of which was determined by the membrane 
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potential. Each region of potassium permeability in the membrane con- 
tained four sites. If all the four sites in such a region were simul- 
taneously occupied by a charged particle, potassium was allowed to move 
through the membrane at that region, otherwise the membrane remained 
impermeable to potassium. In terms of the mathematical formulation, n , 
potassium activation, represents the probability that a charged particle 
is at a site. The rate constant O. denotes the rate of movement of 

the particles onto the sites while BS represents the rate of particle 


movement away from the sites. 


Sodium permeability was also assumed to be determined by 
sodium specific charged particles in the membrane. However, in contrast 
to the situation hypothesized for potassium conductance, there were two 
distinct types of particles, an activating particle and an inactivating 
particle, which acted independently of each other. The membrane became 
permeable to sodium in a certain region when three of the four sites 
contained in that region were simultaneously occupied by activating 
particles. The region became impermeable to sodium if at least one of 
the sites was occupied by an inactivating particle. In terms of the 
mathematical presentation, m , sodium activation, represents the 
probability an activating particle is at a site, and h , sodium in- 
activation, denotes the probability an inactivating particle is not 
occupying a site. O, and BL represent the rate of movement of the 
activating and inactivating particles, respectively, onto a site. B 
and a denote the rate of movement of the respective particles in the 
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From this model proposed for the mechanism of permeability 
changes it is evident the objective of trying to conform to Goldman's 
equation was an attempt to account for the variation of the rate constants 
ay membrane potential in terms of the effect of the electric field on 


the movement of the charged particles. 


The Hodgkin-Huxley equations do not account for the passive 
bahavior of the membrane to electrical stimulation. In the form (1.2.4) 
- (1.2.12), they represent a squid nerve axon with the experimentally 
induced condition of a uniform membrane potential over the entire axon 
length or a "space-clamped" axon. With this assumption, the total 
membrane current is zero except during stimulation. The equations have 
been extensively examined to test their ability to simulate the actions 
of the squid giant nerve axon under these conditions. Although the 
existence of a unique solution has been demonstrated (Cole, Antosiewicz, 
and Rabinowitz: 1955), a closed solution has not been obtained. Thus 
the behavior of the equations has been studied by numerical methods and 


computer simulations. 


The original investigation of the validity of the equations 
was carried out by Hodgkin and Huxley (1952d). To determine the res- 
ponse of the membrane potential to an impulse stimulation administered 
at time t = 0 , Hodgkin and Huxley solved the equations using the 
numerical method of Hartree (1932-3). The form, amplitude, and thres- 
hold of the calculated membrane ees potential compared well with the 


experiment results. By linearizing the equations, the damped oscillatory 
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behavior of the membrane potential resulting from small amplitude 
constant current stimulation was reproduced. The equations also 
successfully predicted the form of subthreshold responses, anode break 
responses, the response to stimulation during a refractory period, and 


impedance changes during the action potential. 


Under certain circumstances the predictions of the model are 
not in keeping with known electrophysiological results. The equations 
predict an indefinite train of impulses for stimulation by a rectangular 
wave current of strength above a critical level (Cole, Antosiewicz and 
Rabinowitz: 1955; FitzHugh and Antosiewicz: 1959). However, for a 
"space-clamped" squid nerve, stimulation by a rectangular wave current 
always results in at most a finite number of impulses (Hagiwara and 
Oomura: 1958). This represents a serious discrepancy between the 
predictions of the model and the behavior of the squid nerve. The 
model also does not satisfy the all-or-none property of the nerve. 
FitzHugh and Antosiewicz (1959) using a digital computer demonstrated 
the existence of a continuous gradation in response amplitude ranging 
from that of a subthreshold response to that of a full action potential. 
This particular difference between the computed and experimental results 
is not considered to be serious but rather the result of the digital 
computer not being programmed to account for the spontaneous fluctuations 


in the membrane potential (FitzHugh: 1969). 


Further testing of the equation's ability to simulate the 


behavior of the squid nerve has required modification of the equations 
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in a manner prescribed by experimental results. Certain experiments 
indicate that the dependent variables of the Hodgkin-Huxley equations 
can be altered by variables not accounted for by the equations. By 
modifying the equations in the manner determined by the experiment, 
the predictions of the equations can be compared with the behavior of 


the nerve. 


Frankenhaeuser and Hodgkin (1957) observed that the membrane 
potential varied with the extracellular calcium ion concentration. The 
relationship between the change in membrane potential AV and the 


calcium concentration [Ca] was experimentally determined to be 


tie [Ca] 
Ve= O.32 in TCa]_ 


where [Ca] | is the normal calcium concentration. Huxley (1959) 
incorporated this in equations (1.2.4) - (1.2.12) by substituting 
V+AV for V.. The predictions of these modified equations agreed 


well with the experimental results for different calcium concentrations. 


By injecting tetraethylammonium chloride into the squid nerve 
axon, an above threshold stimulation produces a long lasting action 
potential (Tasaki and Hagiwara: 1957). FitzHugh (1960) was able to 
reproduce this phenomena with the equations by reducing O., and B. 
by a factor of 100 and increasing O, and By byy avtiaetor. of WiSl. 
These modifications are supported in part by the voltage-clamp data 
which revealed that during the first few milliseconds of the response of 


a nerve injected with tetraethylammonium, the potassium current is 
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either absent or significantly delayed. 


Calculations using these modifications agree with other 
experimental results also obtained from a nerve treated with tetra- 
ethylammonium. These include the abolition of the response to a 
repolarizing current of a critical threshold value applied during the 
plateau phase of the response, and the response to a second stimulation 
following the long lasting action potential. The impedance changes 


occurring during the plateau, however, are not reproduced by the equations. 


Another validation of the equations using an analogous 
approach has been carried out to predict the effect of high extracellular 
potassium concentration on the squid nerve axon (George and Johnson: 1961). 
In this case only partial success was obtained in reproducing the 


behavior of the nerve. 


The Hodgkin-Huxley theory, with certain exceptions, has 
demonstrated a remarkable ability to simulate the behavior of the 
"space clamped" squid axon to electrical stimulation. However, Hodgkin 
and Huxley (1952d) expressed two important reservations regarding the 
formulation. First, the equations have been empirically determined 
and as such may not be the only satisfactory formulation. Second, the 
success of these equations does not necessarily provide a proof that 
the hypothesized mechanism for membrane permeability changes is valid. 
It is in the light of these two remarks that alternative formulations 
and interpretations of the Hodgkin-Huxley theory have arisen. These 


alternative formulations will be discussed in 81.3. 
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81.3 Alternative Formulations for the Hodgkin-Huxley Theory 


The reservations regarding the Hodgkin-Huxley formulation 
expressed in the last paragraph of §1.2 have fostered two avenues of 
research. Investigations have been initiated to determine alternative 
chemical and physical mechanisms for membrane permeability changes which 
satisfy the Hodgkin-Huxley equations, and alternative empirical formula- 
tion which satisfy the experimental data obtained by Hodgkin and Huxley 
(1952d). Although the success of the equations does not establish 
conclusively the validity of the hypothesized mechanism of permeability 
changes, it does lend support. Other physical and chemical mechanisms 
which satisfy the equations therefore must also be considered in the 
realm of possibilities. These hypothetical mechanisms serve to provide 
a guide line for future experimental work to determine the actual 
mechanism. Alternative mathematical models which satisfy the experi- 
mental results could also suggest different physical or chemical 
mechanisms for the membrane permeability changes. In these cases the 
mathematical formulations are available to distinguish between the 
mechanisms. With the mathematical formulations it is possible to predict 
the behavior of the squid nerve axon thereby providing a quantitative 
means to distinguish the models and establish to a certain degree the 
validity of each permeability mechanism. In 81.3, certain of the 
alternative empirical Ep eri geisne dace mechanisms for the permeability 


changes which have been derived, will be described. 
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In 1963, Agin derived an alternative theoretical basis for the 
conductance equations (1.2.5) - (1.2.12) of the Hodgkin-Huxley model, 
thereby providing an alternative physical framework on which to formulate 


the mechanism for membrane permeability. 


The motivation for the model was the form of the equation for 


the rate constant ome resulting from a transformation. By defining 
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In this form O is proportional to the mean value function of 
statistical mechanics. Further developments result from considering the 
solution of the potassium conductance equation (1.2.5) for a given 


constant membrane potential and initial condition 


0 - 1 
(255) n= ae a +e 
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By writing BO = be , where b= = and employing (1.3.2), the initial 


condition (I.3.3) can be reformulated as 
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Int fornseCies.5) and (1.3.7), n, and n. represent statistical 


distribution functions with properties identical to equilibrium distri- 
bution functions of statistical mechanics. Similar formulations can be 
obtained for the rate constants On and BF , and initial and steady 
state values of m. Furthermore the rate constant BL of sodium 


inactivation given by (1.2.12) is analogous to the Fermi-Dirac distri- 


bution function. 


This analysis suggested that a possible theoretical framework 


on which to base the empirical conductance equations could be obtained 
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using statistical mechanics. To illustrate this hypothesis, Agin, 

using an abstract model satisfying appropriate probabilistic constraints 
derived an equation which is formally equivalent to the potassium 
conductance equation (1.3.4). Although the analysis was not carried 
out, it was suggested that in an analogous manner a model could also 

be derived which would produce equations formally equivalent to the 


sodium conductance equations. 


The analysis of Agin suggests a possible framework for a 
physical model of membrane permeability, the mathematical description 
of which is the Hodgkin-Huxley equations. The original Hodgkin-Huxley 
formulation remained unchanged. An alternative form of the potassium 
conductance equation (1.2.5) was derived by Cole and Moore (1960). It 
was found necessary to modify these equations in order to account for 
a time delay in the potassium current I » resulting from the 


depolarization of a prehyperpolarized squid axon. The required 


reformulation is 


(153.8) Be = Sy 0 


A 


(1.3.9) a =o (1-1) -6n 


where &x is a constant, and both O, and B are functions of the 
membrane potential. With the appropriate values. of OL,» Bn and initial 
condition for (1.3.9), the reformulation proved adequate in representing 
Bx for the experimental conditions originally considered by Hodgkin 
and Huxley (Cole and Moore: 1960). The model thereby extends the range 
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of membrane potentials over which the equations accurately simulated the 


behavior of the squid axon. 


A satisfactory alternative formulation for both the potassium 
and sodium conductance equations of the Hodgkin-Huxley model was proposed 


by Hoyt (1963). 


In the presentation by Hoyt, potassium conductance was 


Bx 


assumed to be a function of the variable Ve» which obeys the first 


order differential equation 


dv, 
Cbs3 21:0) ae = Ay Vy. - Vy) 
where Vig oo and ay, are functions of the membrane potential. The 
functions 
Be = By (vy) 
en Vezoo6V) 
and a, = a, (V) 


where V , as in the Hodgkin-Huxley formulation denotes the difference 
between the membrane potential and resting potential, were empirically 
determined in order that for a given membrane potential, the graph of 
Be = g(t) agreed with the experimental results obtained by Hodgkin 
and Huxley (1952d). 


Sodium conductance Bx, was defined to be a function of the 


single variable Va which satisfies the second order differential 
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equation 


Chis ci) 


where Viaes » Y and 6 are functions of the membrane potential. As in 


the case of potassium conductance, the functions 


BNa cf By Na? 

Vac i vee? 
aa) 
6 = 6(V) 

and the initial condition 

d 

Vg 6b »V) a 
dt t=0 


were chosen to obtain the best fit of the curve 8na = Bua ft), for a 


given membrane potential, to the results of Hodgkin and Huxley (1952d). 


In order to simplify the analysis, Hoyt proposed that equation 


(1.3.11) results from the system of coupled first order differential 


equations 

du 

Na 
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where k, for i= 1,2,3,4 are functions of the membrane potential. 


This system of equations is equivalent to the second order differential 


equation 
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For a constant membrane potential, equation (1.3.13) reduces to the form 
of (1.3.11) and the functions k = k, () for. =") 23.4 ‘can be 


determined by comparing coefficients. 


The validity of Hoyt's formulation was examined by computing 
the form of a membrane action potential. Comparison of the result with 
that computed by Hodgkin and Huxley (1952d) revealed a difference only 


during the repolarization phase. 


Hoyt also proposed a possible physical mechanism for the 
membrane permeability changes, which would satisfy the mathematical 
formulation. In the case of Arenas che membrane was assumed to possess 
potassium specific pores,each of which contains a series of p dipoles. 
Each dipole possesses two stable orientations, of which only one allows 
the potassium ion to pass through the membrane. A quantitative 
description of the mechanism for permeability changes was obtained by 
assuming the two stable orientations to be separated by an energy 
barrier e«. ce would therefore correspond to the energy required to 
orient a dipole to permit the passage of a potassium ion,and a pore 
with j4(j = 0,1,2,...,p) dipoles oriented to permit the passage of 
potassium could be considered to possess the amount of energy je . 


With this formulation a pore resembles a quantum system with a ground 
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state, corresponding to the case in which all dipoles are oriented to 
prevent the passage of potassium, and p excited states each determined 


by the number of dipoles oriented to allow the passage of potassium ions. 


A collection of such pores obey Boltzmann statistics in which 
case the probability En ,» that a pore has all dipoles in a potassium 


permeable orientation and the pore is open, is given by 
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where k is Boltzmann's constant and T is the absolute temperature. 
Further, the probability y that in equilibrium a dipole is oriented 


to allow potassium to pass is 


(I. 3.16) i cae ee 


Substitution of (1.3.16) into (1.3.15) yields 
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By comparing the graph f = = 

y comp g graphs of 2. = g,(v,) and De Baye : 
Hoyt determined potassium conductance &, asa function of the 
probability Baty) that a pore is open. Agreement between the two 
functions was attained with p = 9 and the linear transformation 


y = 0.357 + 0.00394v,, . The formulation for defined as a function 


Chen Le) Se = BNP, Cy) 


where g, denotes the conductance of a pore with all dipoles oriented 
to permit the passage of potassium and N is the number of pores per 
unit area of membrane. This formulation enables the determination of the 


range of v thereby providing the means to approximate the delay of 


K b 
the potassium current I. » resulting when a depolarizing voltage-clamp 


is preceeded by a large hyperpolarization. The calculated delay was in 


close agreement with that determined by Cole and Moore (1960). 


Hoyt's hypothesized physical model for the potassium per- 
meability mechanism presents one complication. The difficulty lies in 
the fact that for large changes of Vy the calculated energy barrier 


€ between the two orientations is too large to permit the proposed 


serial arrangement of the dipoles. 


The sodium permeability changes, Hoyt suggested,could be 
modeled in a similar manner except in this case the dipoles would 
possess three stable orientations only one of which would allow a 


sodium ion to pass. 
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A feature common to the models proposed by Hodgkin and Huxley, 
and by Hoyt is that potassium conductance at each permeability region or 
pore is governed by a number of identical events which occur independently 
oF each other. If the potassium conductance equations derived by Cole 
and Moore are interpreted in the same manner as those of Hodgkin and 
Huxley, then twenty-five identical events must occur at each potassium 
specific en the in order for the membrane to become permeable. Tille 
(1965) suggested that it is difficult to reconcile the fact that these 
events occur independently of each other and at the same time are 
localized within a region or near a pore in such a manner that the 


action of any one can prevent the passage of a potassium ion. 


Johnson and Tille (unpublished) determined that the experi- 
mental data of Hodgkin and Huxley (1952d) for potassium conductance con- 
sidered as a function of time for a given membrane potential can be 


empirically described by 


(1.3.19) Be = & 0 
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where k_ and se are constants, and both w and §8 are functions of 
the membrane potential. Equation (1.3.20) lends itself to the physical 
interpretation of describing a system of dipoles each of which can exist 
in one of two stable orientations but which can change orientation only 


if an adjacent dipole is in a specific orientation. 
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This analysis of Johnson and Tille laid the foundation for a 
mathematical model of potassium conductance, proposed by Tille (1965), 
which would avoid the above mentioned difficulty of the previously 
described models. The model avoids the complication by requiring only 
a single event, the reorientation of a dipole or the movement of a 
blocking particle to open or close a potassium specific pore. One 
restriction to this simplified model was however required; that is, a 


pore could change state only if an adjacent pore is open. 


In the mathematical formulation of the model, potassium 


conductance was defined to be 


Bx 


(15.21) Se = Be N 


K 


where By is a constant and n denotes the probability a given pore 
is open. By assuming the pores to be distributed over the axon in a 


Square array then 


2 3 
L¥. 3222) 1 - an(l - n)(4 - 6n + 4n? - n°) - Bn(4 - 6n + 4n - 7) 


where n(1 - n)(4 - 6n + an® ~ n°) is the probability a pore is 

closed and at least one adjacent pore is open, and n-th - 67 + ic ~ oe 
is the probability a pore is open and at least one adjacent pore is open. 
a and 8 denote the rates of departure and arrival of the blocking 
particles, respectively. These rate coefficients were assumed to be 


functions of the membrane potential and were empirically determined to 


be 
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Rearrangement of (1.3.22) reveals that it is of a form closely 


resembling (1.3.20). 


Computations carried out using (1.3.21), (1.3.22) and (1.3.23) 
for the potassium conductance equations in the Hodgkin-Huxley formulation 
gave good approximations of the experimental results obtained by Hodgkin 


and Huxley (1952d) and by Cole and Moore (1960). 


Another model for the mechanism of potassium permeability 
changes which avoids the difficulty of prescribing the precise number 
of events occurring at a region or pore was described by FitzHugh (1965). 
This model assumed that associated with each pore is an infinite number 
of sites. If at least one of the sites is occupied by a potassium 
specific charged particle, the potassium ion is prevented from moving 


through the membrane. 


Potassium conductance By » as in Tille's formulation, was 
defined to be 
Ci3.24) Be = 8 N 


where oe is a constant and n is the probability a pore is open. To 


obtain the dynamic behavior of the system, it was supposed the charged 
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particles enter and leave the sites at random, independently of the 
number of sites occupied. Further, the rate of arrival and departure 

of the particles is governed by the rate constants a and 8 ; 
respectively, both of which were assumed to be functions of the membrane 
potential. With these assumptions, the probability = that a given 


pore has j sites occupied is 


dx 
(PSt25) Se on ee oe De 


j 
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The steady state solution of (1.3.25) for a given membrane potential is 


oJ 
(3) x, 
(173.26) x, =\———=—aform j = 0,1,2,... . 
j j! 
Since ) x, = 1 then n ae ,» the probability a pore is open, is 
(eae as 
e B and equation (1.3.26) can be explicitly written as 
4 mae 
Ge ® 
Chon aJd) x enter jr it Fore ques 0 2 ee. pis 


By (1.3.27), for a given membrane potential the number of occupied sites 
for a given pore in a steady state obeys a Poisson distribution with mean 


a Assuming the number of occupied sites for a given pore also obeys 
a Poisson distribution during the transient state, that is assuming the 
Jer 
wn _ (ure 
solution for a step current applied at t = 0 to be x, = af aque 


transforms (1.3.25) to 
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Therefore, since the derivation of (1.3.28) is independent of the 
functional dependence of a and 8 on membrane potential and time, 


the dynamic behavior of n is given by 


1013729) neelt 


where Uw satisfies (1.3.28). To determine the functions oa = a(V) 
and § = 8(V) , FitzHugh assumed for convenience that these obeyed 
the Eyring theory of reaction rates. The functional dependence of a 


and 8 one the membrane potential was therefore given by 


H 


( ¢c + ari 
| 
| a(V) =Ae = 

G13 330) 
a 
{ ee 
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A, B , and C are constants, H denotes the energy barrier which the 
particles must cross to enter or leave a site, r is the fraction of 
the total membrane potential between the two wells, and N is the 


number of electronic charges on the particle. 


FitzHugh's formulation for potassium conductance is therefore 
given by (113.24). 1.3.28), (1.3229)... ‘and €1.3530).) Withsempirically 
determined values of a and §8 and the appropriate initial condition 


for (1.3.28), the solution of the conductance equations fitted the 
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experimental results of both Hodgkin and Huxley (1952d), and Cole and 


Moore (1960) with reasonable accuracy. 


The basic model can also be generalized to include the case 
in which more than one species of charged particles controls the 
permeability of the ion, as is the case hypothesized for sodium 
permeability in the Hodgkin-Huxley model. For the general model, the 
pore was assumed to be closed if at least one particle, irrespective 
of which type, occupied a site. The formulation was derived by 


moditiying? ICI .36:29)4 and* (Ls 3:130)nuto 
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respectively, where p is the number of distinct species of charged 
particles, and Ol, and B, represent the rate of movements of the 


2 th e 
particles of the i species. 


An interesting feature of the model proposed by FitzHugh is 
that by assuming the number of sites associated with each pore to be 
finite reduces the formulation for potassium permeability to that of 
the Hodgkin and Huxley model or Cole and Moore model. This results 
from the fact that with a finite number of sites, the Poisson distri- 
bution is replaced by a binomial distribution in which case (1.3.29) 


and (1.3.28) would be replaced by 
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where S denotes the number of sites at each pore. For (1 - = =n, 

S = 4 , and the appropriate definition of a and 8 , equations (1.3.24) 
(1.3.33), and (1.3.34) give the original formulation for potassium 
conductance proposed by Hodgkin and Huxley. With S = 25 the equations 


are formally equivalent to the model of Cole and Moore. 


A number of different models have been discussed in 81.3 and 
in each case certain aspects of the behavior of the squid nerve axon is 
accurately simulated by the models. There do, however, exist certain 
similarities between the models. In particular, the mathematical 
structure of the equations for certain of the models are similar in 
their formulation for certain qualitative features of the axon necessary 
for excitation. These structural similarities permit the classification 
and comparison of the models. Before considering such a classification 
of the models it is desirable to extend the variety of excitation 
models considered. A class of models which are formulated in a different 


manner from those discussed to this point will be described in $81.4. 
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81.4 Physical and Chemical Models for Nerve Excitation 


For the mathematical models of nerve excitation considered in 
81.2 and $1.3, the equations were empirically determined then a physical 
mechanism, which satisfies the mathematical formulation, was porposed to 
explain the membrane permeability changes. A number of mathematical 
models for excitation have been derived by postulating the mechanism of 
permeability changes the then deriving the equations. In each case, the 
meahcnism was based on certain physiological or histological features of 
the squid axon which are conceivably related to such permeability changes. 
The mathematics was then employed to describe the model and make predic- 
tions regarding it. In 81.4 three models which were constructed in this 


manner will be considered. 


Mullins (1959) proposed that the ionic currents result from 
the movement of ions through membrane pores. In contrast to previously 
discussed pore models, it was supposed that the pores are not ion 
specific. Ionic selectivity of the membrane was accounted for by 
assuming membrane pores discriminate between ions on the basis of ion 
size. For the model, a permeable ion is permitted to pass through a 


pore only if the ion and pore have identical radii. 


The membrane permeability changes were attributed to changes 
in the distribution of the pore sizes. At zero membrane potential, the 
pore sizes were hypothesized to be distributed with a Gaussian distri- 


bution such that the mode of the distribution corresponds to the radius 
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of the potassium ion. Reduction of the membrane potential to that of 
the resting membrane potential results in the continuous transformation 
of the distribution of pore sizes present at the zero membrane potential 
to a Gaussian distribution with mode that of the radius of a sodium ion 
at the resting membrane potential. This change of the distribution of 
pore sizes with membrane potential was assumed to follow from the 
deformation of the membrane by nonpermeable ions and ions of low 
permeability. Such ions are forced into the entrances of certain pores 
by the non-zero membrane potential, thereby altering the size of the 
pore to that of the intruding ions. The ions considered to be responsible 
for this action are calcium ions at the extracellular membrane surface 
and isethionate ions at the intracellular surface. In bouh cases the 
ionic radii are equivalent to that of the sodium ion and at the resting 
membrane potential these ions are assumed to block all pores of that 
radius. Subsequent depolarization of the membrane removes the intruding 


ions from a blocked pore and restores the pore to its undistorted size. 


The axon membrane under the influence of an electrical field 
exhibits the properties of a parallel plate condenser. In this case the 
mechanical force exerted by the calcium or isethionate ion is proportional 
to the square of the membrane potential. Since this applied force results 
in the redistribution of pore sizes, Mullins proposed that the mode r 
of the distribution of pore sizes is a linear function of the square of 


the membrane potential given by 
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where Ty and Tyq are the radii of the potassium and sodium ion, 
respectively, and k is a constant. With (1.4.1), the probability ho 


that the extracellular entrance to a pore is of the sodium ion radius 


at a given membrane potential is given by the error function 


reryté 


(e-r)? 
(1.4.2) h =f ces "int dx 


where 6 is an arbitrarily chosen tolerance and c is aconstant. A 
similar formulation can be written for h, ,» the probability the intra- 
cellular entrance of a pore is of the sodium ion radius. By assuming 
that in the steady state all sodium sized pores are blocked at both 
entrances, Mullins succeeded in deriving expressions for both potassium 


and sodium conductance in terms of hy and h, $ 


For this physical model proposed for ion permeability and 


membrane permeability changes, potassium conductance By is defined to 


be 
CT 24,3) Be = Sy I,d4 


where ip and are the probabilities that extracellular and intra- 


Jy 
cellular entrances, respectively, of a pore have the radii of the 
potassium ion, and ees is a constant. Calculations indicated that the 
results of this formulation were comparable if Hi and 4, are replaced 
by b+ (1-h.) and b+ (1-h,) , where b denotes the probability that 


the entrance of a pore is of the potassium ion size in a fully distorted 


membrane. The formulation for potassium conductance was therefore given 
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by 


(1.4.4) Bere TN did nc 


The expression for sodium conductance is more complex 


eNa 
than that for potassium since it is necessary to account for the compe- 
tition between sodium and the blocking ions, and the inactivation of the 
sodium current. In the model, the blocking ions were assumed to be 
embedded in the pore entrance in a partially hydrated form. Removal of 
such an ion from the pore requires energy to completely hydrate the 
ion. The effect of this energy barrier at the extracellular surfaces 

was incorporated in the formulation by supposing that upon instantaneous 
depolarization, the distribution of entrances which remain blocked is 

an exponential distribution with mean ko . At the intracellular surface 
the distribution is exponential with mean k, . Inactivation of the 
sodium current was assumed to be the property of the membrane itself. 
Once the intruding ion is removed from the entrance, the internal 
structure of the membrane restores the pore to its undistorted size. 
Hence sodium permeability is prevented if the undistorted size of the 
unblocked pore is not that of the sodium ion. Mathematically sodium 
inactivation was expressed in terms of the distribution of pores which 
have been unblocked by depolarization and have returned to their un- 


distorted size. The distribution was chosen to be an exponential with 
mean k ‘ 


Thus for sodium permeability to occur, a pore of sodium ion 


radius must become unblocked and a sodium ion must penetrate the pore 
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before it attains a size outside the limits which permit the passage of 
this ion. The equation chosen to describe the probability the entrance 
to a pore is both unblocked and of the sodium ion size was that for the 
concentration of an intermediate in a series of chemical reactions. 


Sodium conductance By, was therefore given by 


Cl. A...) 


= Ah 
0) 


®Na 8Na ; 


rae 
} 
} 
: 


i 


7 


-k,t “kt | 
e (3) 


=? 


where ras is a constant, t denotes time, and Ah , and Ah, represent 
the probabilities an extracellular and intracellular entrance, respectively, 
is open during the depolarization. The quantities in square brackets of 
(1.4.5) represent the probabilities an entrance is both unblocked and of 


the sodium ion size. 


y 


The mathematical description for the electrical behavior of the 
squid axon proposed by Mullins consists of (1.2.3) - (1.2.5) of the 
Hodgkin-Huxley equations together with (1.4.1) - (1.4.5). Although 
extensive analysis of this formulation has not been carried out it has 
been demonstrated that the formulation for potassium conductance consi-_ 
dered as a function of membrane potential yields accurate results 


(Mullins, 1959). 


For Mullin's theory, the structure of the membrane was 
significant only in the inactivation of the sodium current. To incorporate 


this membrane phenomena in the theory it was not necessary to provide a 
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description of the membrane structure. A model for nerve excitation 
which requires the explicit description of the structure of the membrane 
was proposed by Goldman (1964). The model was based on the electron 
microscope evidence which indicated the axon membrane consists of a 
bimolecular layer of lipid covered on either side by a protein layer. 
At intervals throughout the bilayer, phospholipids are also found to be 
present. Each phospholipid consists of a polar phosphate end group, 
which is situated at the surface of the membrane, and a hydrocarbon 
chain, which extends into the interior of the membrane. A polar end 
group which is not bonded to a membrane protein forms a dipole which 
can bind to ions in the external medium. This fact forms the corner- 
stone of Goldman's theory for cationic permeability. In order for a 
cation to penetrate the membrane it must initially be bonded to such 

a dipole. Permeability is initiated by the collision of an unbonded 
cation with the bonded cation with sufficient energy to break the bond 
and drive the released cation into the membrane. The released cation 
then diffuses through the membrane interior which is assumed to be a 


continuous, inert medium. 


To account for the ionic specificity of the membrane, Goldman 
proposed that an unbonded phosphate group, which can change configuration 
with the alteration of the electric field and ionic environment present, 
would in this manner change its combining properties for different ions. 
For each phospholipid, three interchangeable ion specific configurations 
of the phosphate end group were considered. These are denoted by Cys 


Under the influence of the electric field for the 
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resting membrane potential, configuration Cy » which favors the bonding 
of calcium ions, predominates. Removal of this electric field by 
depolarization results in a predominance of dipoles in configuration 

Cir » for which the phosphate group has an affinite for sodium ions. 
Also present in a significant proportion during depolarization of the 


membrane is the potassium specific configuration C By considering 


i ae 


the interchange between Cir and Clit to be independent of the electric 
field present, Goldman proposed that during depolarization sodium 
permeability would attain a maximum then decline and be replaced by 


potassium permeability. 


If it is assumed the ionic exchange at the polar end group 
occurs with sufficient rapidity that the configurational changes repre- 
sent the rate-limiting steps it is possible to quantitatively describe 
the ionic currents resulting at the membrane surface in terms of the 
end groups and configurational changes. For convenience in determining 
the mathematical formulation only phosphate end groups at the extracellular 
surface membrane were assumed to participate in the membrane transfer of 
cations. To obtain the quantitative description, Goldman considered the 
total number n,, of end groups per unit area of extracellular membrane 
surface which are in the form of dipoles. Of this total My» Myy > and 
Mr denote the number with configuration Ci» Crt » and Crit ; 
respectively. The ionic currents were defined in terms of these variables 


and their dynamic behavior which is dependent upon the interchange of the 


configurations and the effect of the ionic environment present. 
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The rates of interchange among the configurations are governed 


by the rate constants k., k., k 


ie eZ 


k d d 
1 an k. enote the rate of change from C, into Cir and Chir 5 


respectively. k, and ke represent the respective reverse transfor- 


mation rates. k denotes the rate of change of end groups from 


configuration Cir into Crit » while ky denotes the rate of change 


in the opposite direction. Since ki» ky» k, and ke are dependent on 


the electric field strength, it is possible to derive explicit expressions 


3? Kye k. » and ke . Rate constants 


for these in terms of the energy required to change the configuration. 


The formulations are 


yO 
f = s 
| ky boa bi 
oil 
k, = mae 
livdoes 2 
(124.6) 4 ae 
| k_ = me 
i 5 9 
ay, 
: k,. = moe % 
6 6 


v6 -y6 
where e ~ and e ® are expressions which account for the energy 


required to change the configuration, of is the potential across the 


surface of the membrane, and Ms My, M, and m, are constants. 


The effect of the ionic environment present is to restrict the 
number of phosphate end groups which are available for configurational 
changes. In particular, end groups with configuration Cy or Crit 
which are bonded to their preferred cations can not participate in 
configurational changes. In order to determine the number of end groups 


which are bonded to cations, it is necessary to take into account the 
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velocity and concentration of unbonded cations and the number of 
available end groups. For those phosphate end groups which are in 
configuration Cy » the number Mo, per unit area of membrane which are 


bonded to calcium is given by 


= v 
(i Ae) Mo, = P [Ca] at 


where [Ca] , is the concentration of calcium in the extracellular 
medium, p is a constant representing the velocity of ions, and ny 
is the number of end groups with configuration Cy which are not 
bonded to a calcium ion. The corresponding formulation for end groups 


in configuration C is more complex since potassium ions are present 


III 
in the interior of the membrane as well as the extracellular mediun, 
which is not the case for calcium. The number ny of end groups per 


unit area of membrane with configuration Crit which are bonded to a 


potassium ion is 


Cae) a (p (Kl, + p, [KI], 2544 


where [K] and [K], are the concentrations of potassium in the 
extracellular medium and membrane interior, respectively, Be and Py 
denote the velocities of the ions in the respective media, and Mr 

is the number of unbonded end groups with configuration Coir . Those 

end groups with configuration Cir were assumed to be weakly bonded 

to the sodium ions and therefore are not affected by the ionic environment 


present. 


Having established a mathematical description for the configu- 
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rational changes and the effect of the ionic environment, the dynamic 


behavior of the ion exchange process occurring at the end groups is 


given by 
dn, 
—SP = = t q 
| dt (ky + kp)n, + kony, + ken y 
} dn 
(1.4.9) C ee 
de Sehr TE gtir -— t+ ke) tyre 


= = pee v v 
Nyy = Oya ~ Bp 7 (Oy + aK, + By + Ry) 
together with (1.4.6) - (1.4.8). 


The ionic currents resulting from the transfer of ions at the 
end groups was determined in terms of the inward and outward ionic fluxes. 


In the case of potassium, the inward ionic flux is given by 


_ “Ke 
> kT 
Cr, 4.410) ry = on [K] , e 
and the outward flux is 
abe St 
+ kT 
(Dhl) Wee Orn (KI, e 


where Wy and Wey represent the energies required to displace a 
e 

bonded potassium ion, k is the Boltzmann's constant, T is the 

absolute temperature, and both Ove: and On, are constants which 


account for molecular velocity and membrane area factors. The resulting 


formulation for the potassium current I. is 
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Eee Ke | 
i kT i Th 
ees 2 = { aaa 5 
( 12) I, qx yy IKI e Ox l[K] é | 


where q is the electronic charge. Similarly the sodium current 
resulting from transfer of ions at the extracellular membrane surface 


is given by 


Che4el3) Ia = qq | Onait Nal se 


Since the membrane interior was assumed to be inert and 
continuous, upon reaching the interior the ion is only under the influence 
of electrical and concentration gradients. To this point only the 
permeability of cations has been considered. In the case of anions, 
permeability was assumed to be independent of the phosphate end groups 
and resulted entirely from diffusion. The permeable anion is therefore 
only affected by the electric and concentration gradients across the 
membrane. Thus the formulations for the ionic current of the anions, 
predominately chloride, and that for the gations in the membrane interior 


are analogous and given by the flow equations 
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where 9 is the potential across the interior of the membrane, and Qs 
Oa » and ond are the diffusion coefficients of the respective ions 


in the membrane. In addition, the continuity equation and Poisson's 


equation are also satisfied by the system. Hence 


ox q ot 
OL 
Na _ d[Na] 
CL. #eb5) ope Pama at 
Oley venga iCLi 
Ox 4 dt 
and 
y) 2 
(1.4.16) i= af on ([K] + [Na] - [C1]) 
Ox 


where D is the dielectric constant for the membrane interior. 


In order to examine the formulation presented by Goldman 
it is necessary to obtain approximations for the energy terms. Although 
the analysis is limited, approximate solutions show good agreement with 
experimental data for steady-state current-voltage relations and potassium 


tracer flux ratios. 


The final model to be discussed in 81.4 was proposed by 
McIlroy (1970). In this theory the membrane permeability changes are 
determined by enzymatic activity. The model is of considerable interest 


since the formulation proposed for potassium conductance is formally 
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equivalent to that proposed by Hodgkin and Huxley (1952d), even though 
the mechanisms for membrane permeability changes in the two theories 


differ considerably. 


McIlroy's model represents the extension of two 
previously obtained results. Bass and Moore (1968) 
observed that a critical depolarization of the squid axon membrane 
produces an increment in the membrane pH . Furthermore, an increment 
in membrane pH of such magnitude was found to be sufficient to initiate 
an action potential by nonelectrical means. In an attempt to associate 
changes in membrane pH with the changes in membrane permeability, Bass 
and McIlroy (1968) considered the effects of pH changes on enzymes. 
The model adopted for study was that of Michaelis-Davidsohn. For this 
model the enzyme possesses a neutral base and a singly charged base both 
of which can accept at most one proton in the formation of a cationic 
acid and neutral acid, respectively. The arrangements of the enzyme are 
denoted by E where Po and Pn represent the number of protons 
on the cout te eae singly charged bases, respectively. Of the four 
possible arrangements, E10 was assumed to be the only active form of 
the enzyme with respect to the enzymatic substrate. Bass and McIlroy 


succeeded in determining the fraction f of enzymes in active form as 


a function of the membrane pH . The formulation is given by 
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124,18 = ve 
( ) Y=) (pH PH aD 10 
and 
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K. 2 
Gi 419) a=2 a 
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The symbols [E, 9! and LE od denote the concentration of the active 
enzyme form and of all the enzyme forms, respectively, K. and Kh 

are the dissociation constants of the cationic and neutral acid, respec- 
tively, and y is a measure of the ratio of proton concentration to that 
at which f = ate . Using (1.4.17) - (1.4.19), Bass and McIlroy 
demonstrated that in the membrane medium the activity of suitably 

chosen enzymes undergo significant changes following the increment in 
membrane pH capable of producing an action potential. McIlroy 

extended these results by proposing a physical model in which enzymes 


under the influence of the membrane pH produce the membrane permeability 


changes. 


The model for nerve excitation considered by McIlroy closely 
resembles that proposed by Nachmansohn (1959). For this theory, 
acetylcholine S , present in the membrane, acts as the substrate for 
an unidentified receptor enzyme R and the esterase enzyme H. At 
the resting membrane potential, acetylcholine is assumed to be stored 
in a lipid-bound form free from hydrolysis. Upon depolarization, the 
active form Me of the receptor enzyme combines with acetylcholine 


10 


to form a receptor protein complex with the configurational form Oana : 


This configuration is reversibly converted to the configurational form 
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Pleska which produces the permeability changes responsible for the 
rising part of the action potential. Hydrolysis of the configurational 


H 
form C, of the receptor protein by the active form E19 of the 


pen 


esterase results in the inactivation of the ionic currents. The dynamic 
behavior is determined by rate constants which govern each step of the 
procedure. The capture of acetylcholine by the receptor enzyme was 


assumed to proceed at a rate determined by the constant k -Trans- 


1 e 
formation of the configurational form C into C is determined 
open closed 


while the reverse rate is given by Recs ; 


conf conf 
The rate of hydrolysis is governed by the constant k, : 


by the rate constant 


A quantitative description of potassium conductance was. 
obtained by considering the membrane to contain specialized potassium 
specific permeability regions. In each region, potassium ions are 
permitted to pass through the membrane only if the region contains N 


receptor proteins of the configurational form C . Hence potassium 


closed 
conductance By is given by 
1.4.20 = g, [C ie 
Se) Bx ~ 8x '’closed 
where [C ] is the concentration of the closed configuration of 
closed 


the receptor protein, and ee is a constant. By considering the rate 
constants governing the process it is possible to determine the dynamic 


behavior of [C ]. The formulation for the dynamic behavior is 


closed 


mae 
conf dt 


ree J 


(1.4.21) u lCriosed: hs Kone closed: open 
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d rR aiden | 
( ) de roe] = kyES)) Fig] — &) Eqg|(C,.en! 
where 
4,23 = 
( ) Moe! DC emett! in Le teeed: 
PR 
A igs esate: 
(1.4.24) Reo, 8. tot tot 
Lo ee R 
if Siti + “a cosh (‘y) 
. pee oe 
ae LL toty 
Ci 25) EL = pin " ; 
1+ Te a cosh ( y) 
and 
ty = (pH - PH ay) 10 
(CL, 4.26) 1 
‘PE 
af K, 
a=2/>—|, i-=R,H 
iy 
n 
The equation for ian given by (1.4.24) differs from that proposed by 


(1.4.17) since it is necessary to account for the effect of saturation of 
the receptor enzyme. This is accomplished by replacing [ Ete as 


Res et 


proposed in (1.4.17) by oe rola 
Certain physiological results permit the simplification of the 
formulation for potassium conductance. Examination of the time constants 


of voltage-clamp data associated with potassium conductance reveal it is 
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possible to consider oe lt ae O . Hence (1.4.21) is replaced by 
ie 2d) Le nen) "it cont LC ieeed) : 
By letting 
* se en! 
(1.4.28) i) ook 
Fi0! 


equations (1.4.20) and (1.4.22) are respectively transformed to 


NN} iN 
_=—~fR } 1 \ *N 
(1.4.29) Be = By I E10! ki | on 
conf | 
and 
d ip a * kK * 
(1.4.30) Se = oid Clint jn} -8Bn 
t K 
conf 
where 
¥ k,[S] 
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| r R | Lot 
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4 K 
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and 
cH 7 
kj E 
* ZU TtOe 
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McIlroy also determined that Rate >> 1 in which case (1.4.30) is 
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further simplified to 


(1.4.33) Gla el ie 


With N=4 , (1.4.29) and (1.4.33) are equivalent to (1.2.5), the 
equation proposed by Hodgkin and Huxley (1951d) for potassium conductance. 
This result is encouraging since no attempt was initially made to simulate 
the formulations obtained by Hodgkin and Huxley in the derivation of this 


enzyme model for nerve excitation. 


The formulation for sodium conductance is identical to that 
for potassium conductance. Hence the same formulation (1.4.29), (1.4.31), 
(1.4.32), and (1.4.33) must be able to describe the different behaviors 
exhibited by the two conductances. The difficulty is resolved by 
assuming the enzymes associated with each conductance have different 
locations in the membrane. Those for potassium conductance are located 
at the interior of the membrane while those for sodium are located at 
the exterior of the membrane. By assuming that the values of the time 
constant for the return of the membrane pH to its resting level 
following a voltage-clamp depolarization differ significantly between 
the two regions, the differences between the behavior of potassium and 


sodium conductances are realized. 


The models discussed in 81.2, 81.3 and 81.4 do not exhaust the 
possibilities which exist for mathematical models of nerve excitation. 
These models do, however, indicate the variety of approaches that are 


used in deriving the models. It has been suggested in §1.3 that certain 
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of the excitation models can be classified in terms of the structure 

of their mathematical formulation. Since it is possible to classify 
the models in terms of their mathematical formulations it may also be 
possible to determine a generalized mathematical formulation for nerve 
excitation. To obtain such a generalization it is necessary to examine 
the mathematical properties of the formulations for the different 
models. Much of the research on the mathematical structure of the 
formulations for nerve excitation models has been done by FitzHugh 


and will be reviewed in 81.5. 
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81.5 Classification and Generalization of the 


Mathematical Models of Nerve Excitation 


The fact the mathematical formulations provide a means to 
establish the validity of the proposed permeability mechanisms of each 
model demonstrates their usefulness in the study of the excitatory behavior 
of nerve. The importance of the mathematical formulations however is 
not restricted to this application alone. As pointed out in §1.3, by 
considering the structure of the equations it is possible to categorize 
the models on the basis of their formulation for certain of the basic 
qualitative features necessary for excitatory behavior. Employing the 
quantitative nature of the formulations in the investigation of these 
categories hopefully will provide a means to determine the exact nature 
of these basic features. In 81.4, this proposal was extended by 
suggesting that through the investigation of the structure and mathe- 
matical properties of the formulations it may be possible to derive a 
generalized formulation for the basic interactions required to describe 
the excitable behavior of the squid nerve axon. 81.5 is devoted to the 
classification, where it is possible, and generalization of the excita- 


tion models. 


It is possible to classify certain of the models into two 
distinct categories on the basis of their formulation for sodium con- 
ductance. In the Hodgkin-Huxley theory, sodium conductance is a 


function of two independent variables which are each dependent on wholly 
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independent processes. For the theories proposed by Mullins (1959), 

Hoyt (1963), Goldman (1964) and McIlroy (1970) sodium conductance is a 
function of one variable which itself is dependent on two or more 

coupled processes. This qualitative mathematical distinction between the 
two categroies has been shown by Hoyt (1968) to lead to quantitative 
differences in their predictions of sodium conductance for different 
voltage-clamp potentials following a conditioning voltage-clamp. Con- 
clusive experimental evidence, however, has not to this point been 
obtained which indicates which of the categories simulates the axon most 


accurately. 


Certain of the models can also be categorized in terms of the 
interdependence, demonstrated in the formulation, between the sodium and 
potassium currents. The models proposed by Hodgkin and Huxley (1952d), 
Hoyt (1963) and McIlroy (1970) consider the sodium and potassium 
permeability changes to be independent of each other. The models of 
Mullins (1959) and Goldman (1964), however, consider the sodium and 
potassium permeability changes to be coupled. This coupling results 
from the assumption that both sodium and potassium ions share a common 
pathway through the membrane. Mullins (1968) conjectures that the 
hypothesis of a common pathway would be invalid if circumstances exist 
for which the sum of sodium and potassium conductance at a given time 
exceeds that of the maximum sodium conductance. Again, experimental 


evidences has not completely resolved which of the proposals is valid. 


In order to obtain a generalized formulation for nerve 
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excitation, it was pointed out in 81.4 that the mathematical structure of 
the existing formulations must be investigated. FitzHugh, in a series of 
papers (FitzHugh, 1960 and 1961) analyzed the Hodgkin-Huxley equations in 
an attempt to determine the physiological importance of each of the 
dependent variables. The analysis was carried out using phase space 


methods. 


To determine the significance of the dependent variables of the 
Hodgkin-Huxley equations and to reduce the four dimensional (V,m,h,n) 
phase space to one which can be graphically portrayed, FigzHugh (1960) 
reduced the system of equations by considering one or more of the 
dependent variables to be constant. It is possible to classify the 
variables V, m, h and n_ on the basis of the magnitude of their 


relaxation times which are given by 
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respectively. Since the magnitude of Th and tT were determined to 


be approximately one-tenth the magnitude of ue and Th» FitzHugh 
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- considered the Hodgkin-—Huxley equations with n and h _ held constant at 


their resting state values. This reduced system of Hodgkin-Huxley 


equations is 


ne 3 oe 
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where no and ho denote the values of n and h at the resting 


membrane potential, V=0. 


The (V,m) phase plane of the reduced system (1.5.2) was 
determined by analogue computer. Trajectories in the plane represent 
the response to various magnitudes of instantaneous current impulse 
administered at t = 0. The reduced system admits only three singular 
points. Two of the singular points were determined to be stable. The 
third singular point, a saddle point, together with two stable separatrices 
form a boundary in the phase plane separating the two stable singular 
points. The strength of the initial current pulse determines on which 
side of these separatrices a trajectory lies and hence to which stable 


singular point it will approach. It is possible to correspond these 
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phase plane regions and points to physiological states of the squid nerve 
axon. The two stable separatrices together with the saddle point 
correspond to the threshold phenomena present in the squid axon, while 
the two stable singular points represent the physiological states of 


rest and excitation, respectively. 


By removing the restrictions on n and h and considering the 
(V,m) phase plane for the complete Hodgkin-Huxley equations produces a 
number of changes. Only the stable singular point corresponding to the 
resting state remains. The threshold phenomena in this case is repre- 
sented by a trajectory, for which only a portion of its length acts as 
a separatrix dividing the neighboring trajectories, which represent the 
responses to instantaneous current impulses, into subthreshold and full 
action potential responses. It was impossible, however, to find a unique 
separatrix which would distinguish between all the trajectories on the 
basis of this property. The (V,m) phase plane, therefore, demonstrates 
the inability of the Hodgkin-Huxley equations to satisfy the '"all-or-none" 
law as demonstrated by the squid nerve axon. In contrast to the (V,m) 
phase plane with h and n constant, the (V,m) phase plane repre- 
senting the complete Hodgkin-Huxley equations contains regions which 
correspond to the physiological states of recovery and refractoriness. 
Thus the pair of dependent variables V and m can be considered to 
represent the property of excitability of the axon while n and h 


represent the properties of recovery and refractoriness. 


The Hodgkin-Huxley equations were considered by FitzHugh to 
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be a member of a class of non-linear systems which demonstrate excitable 
and oscillatory behavior. To further analyze the properties of the 
Hodgkin-Huxley equation, FitzHugh derived a mathematical model which would 
represent this class of non-linear systems and at the same time be more 
tractable than the Hodgkin-Huxley equations. The model was not designed 
to represent experimental data accurately but rather to represent the 
basic interactions required to produce the various physiological states 


of the squid axon. 


The model proposed by FitzHugh is based on the van der Pol 
equation (van der Pol: 1926), for a relaxation oscillator given by 
dx 2 dx 
Mine) —~+c(x - 1) ae x = 0 


Aes 


where c is a positive constant. Applying the Liénard transformation 


(Liénard; 1928) 


Wo 


dx x 
CEs5. 4) y = Bein oa Lae 


to (1.5.3) yields the set of differential equations 


3 
ot = e(y +x - 4) 
| dt 3 
(1.5.5) 
iy * Re 
dt a 


\ 


The phase plane of (1.5.5) contains an unstable singular point at the 
origin. To remove this unstable singular point, FitzHugh considered the 


equations 
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dx x 
dt c(ly +x - Ta7t z) 
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and both a and b constant. In order to consider system (1.5.6) as a 
model for nerve excitation it was assumed that x denotes the displace- 
ment of the membrane potential from the resting potential and z denotes 
the membrane current density. System (1.5.6) is referred to in the 


literature as the Bonhoeffer-van der Pol equations. 


FitzHugh analysized the Bonhoeffer-van der Pol equations by 
phase space methods using computer simulations. Conditions (1.5.7) had 
been chosen in order that for z= 0 _, the isoclines corresponding to 
S = 0 and x 0 intersect at a unique point, which is stable. This 
singular point represents the resting state of the squid nerve axon. It 
was also possible to correspond other physiological states of the axon, 
such as the threshold phenomena and refractory period to specific regions 
of the phase plane. Furthermore, the electrical behavior of the axon 
under a variety of stimulus conditions could be adequately interpreted on 


the phase plane. However, for a step current of magnitude exceeding a 


critical value the trajectories in the phase plane approach a stable limit 
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cycle. Hence, the Bonhoeffer-van der Pol model produces an infinite train 
of impulses for constant current stimulation above a critical level, as 
did the Hodgkin-Huxley model. As pointed out in §1.2, this is contrary to 


the known behavior of the squid nerve axon. 


The reduced system of (1.5.6) with y held constant was also 
considered by FitzHugh (1961). The x phase line of this reduced system 
contains three singular points. Two of the singular points are stable 
and represent respectively, the physiological states of rest and 
excitation. The third singular point, a saddle point situated between 
the two stable singular points, represents the threshold phenomena of the 
axon. There are, however, no regions in the x phase line which 
correspond to the recovery or the refractory period of the axon. The 
dependent variable x , therefore, represents nerve excitability and 
corresponds to the variables m and V of the Hodgkin-Huxley equations, 
while the variable y represents the states of recovery and refractoriness 
and corresponds to the variables n and h of the Hodgkin-Huxley 
equations. This final result provides a means for comparing the 


formulations proposed by Hodgkin and Huxley and by FitzHugh. | 


In order to compare the Hodgkin-Huxley equations and the 
Bonhoeffer-van der Pol equations it is necessary to obtain a phase plane 
of the Hodgkin-Huxley equations which can be compared to the (x,y) 
phase plane. To construct such a phase plane, FitzHugh utilized the 
previously obtained result which ie set the correspondence between the 


dependent variables of the two formulations. The resulting phase plane 
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is the (u,w) phase plane for which u is a function of m and V F 

and w isa function of n and h. The points (u,w) in the plane 
were determined in the following manner. Each value u corresponds to 
the points in the (m,V) plane for which u = V - 36m while each value 
w corresponds to the points in the (n,h) plane for which w = = (n - h) 
The projections in each case had been chosen in order to preserve the 


physiologically significant features represented on the respective planes. 


The (u,w) phase plane of the reduced system of Hodgkin-Huxley 
equations, as with the (x,y) phase plane of the Bonhoeffer-van der Pol 
equations, contained regions and points which eoerettonl to specific 
physiological states of the nerve. Furthermore, the (u,w) phase plane 
and the (x,y) phase plane were found to be qualitatively similar for 
a variety of stimulus current inputs. In particular, for both models 
the trajectories representing the response to stimulation by a step of 
sufficient magnitude, approach a stable limit cycle. These similarities 
Suggest that the Hodgkin-Huxley and Bonhoeffer-van der Pol models can be 
considered as belonging to the same general class of excitable-oscillatory 
systems. FitzHugh (1961), however, expresses one reservation in making 
this correspondence between the two models. The objection is that the 
(u,w) plane is not a phase plane since in general there exist an in- 

“ and seat for each point (u,w) . In 


dt dt 


spite of this reservation, the ability of the Bonhoeffer-van der Pol 


finite number of values of 


model to qualitatively describe the physiological states of the nerve 
and accurately predict certain excitable neural phenomena secures its 


relevance in neural modelling. 
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In 1969 FitzHugh attempted to derive a generalized mathematical 
model of squid nerve excitation. The basis of the mathematical 
formulation is the assumption that the squid axon membrane can be 
considered to consist of a capacitor in parallel with a conductor. Hence 


the basic equation is 


AV taps 
(5. 8) ae” G [I - 1,] 


moere. IT, V, Cu » and t are as defined for the Hodgkin-Huxley equations 
and I, is the current density per unit area of membrane through the 
conductor. Experimental evidence (Hodgkin and Huxley 1952) indicates 

the excitatory behavior of the axon can be attributed primarily to the 
behavior of the conductance current I, . For the generalized formu- 
lation to exhibit the excitatory behavior of the nerve, FitzHugh proposed 
that I, must be a function of variables which, as in the case of the 
dependent variables of the Hodgkin-Huxley and Bonhoeffer-van der Pol 
formulations, represent specific physiological states of the excitation 
process. The functional form of I, was derived from the fact it is 
possible to classify the physiological states or events of excitation 
according to their duration time. The classifications considered by 
FitzHugh, listed in order of increasing time duration, are excitation 
events, recovery events and adaptation events. Excitation events include 
the threshold phenomena, and the rising phase of the action potential, 
while recovery events consist of the falling phase of the action potential 
and the refractory phase. Adaptation events, which develop the slowest, 


account for the reduction of the frequency and possible termination of a 
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train of impulses produced by a prolonged constant stimulation. The 
excitatory nerve behavior was, therefore, incorporated in the formulation 
by assuming I, is a function of the membrane potential and one or more 
time dependent variables each of which can be classified according to 
the system devised for classifying physiological events of excitation. 


Thus the generalized mathematical formulation for nerve excitation is 


given by 
pedy 2 ay 
eure eee haem toe CV We Wo ames Wl 
| dt Fis dtirdoel 242. Walthe Se 
(r-5-9) } 
| ae = F,(VjWy Wy ye+e5W)5 se LSS eee wath, od 


where each be for j =1,2,...,u is either an excitation variable, 


recovery variable, or adaptation variable. 


The mathematical models discussed in sections §I.1 through 

81.5 with the exception of the model of Young (1937) for b #0 can all 
be described by the system (1.5.9). It is of interest to classify, in 
accordance with FitzHugh's proposal, the variables of certain of these 
excitation models discussed. For the models of the two-factor theory of 
excitation, L is a recovery variable while K represents the membrane 
potential. For the model proposed by Hodgkin and Huxley (1952d) m is 
an excitation variable, and both n and h are recovery variables. In 
the case of the Bonhoeffer-van der Pol formulation (1.5.6) which can be 


transformed to the equivalent formulation 
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with a, b , and @ constant and with the conventional notation for 
membrane potential and membrane current in excitation models, W is a 
recovery variable. Neither the Hodgkin-Huxley formulation nor the 
Bonhoeffer-van der Pol formulations contain adaptation variables. 

FitzHugh (1969) suggests that the absence of such variables is responsible 
for the formulation's production of an infinite train of impulses for step 


current stimulation of sufficient magnitude. 


FitzHugh (1969) has attempted to extend the Hodgkin-Huxley and 
Bonhoeffer-van der Pol formulations to include adaptation variables. To 
obtain the required formulation it is necessary to reconsider the two- 
factor theory of excitation. As stated in 8I.1, the two-factor theory 
with only one exception is. unable to account for events following the 
initiation of an action potential. This exception, derived by Katz 
(1936), is the extension of Hill's model (Hill: 1936a) to enable the 
determination of the duration of a train of impulses. Katz proposed 
thatere “Ki <ce) in “(1eS ) dnd AL.1.7 +), then for stimulation by 2 step 
current of sufficient magnitude the local potential K would remain 
above the threshold L for a period of time during which repetitive 
impulses would occur. The effect of the assumption that k << A is to 
change variable L to an adaptation variable. This assumption also 


produces changes in the formulation proposed by Hill (1936a). For 
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k << X equation (1.1.5) is replaced by 


Olek L) K = Ko ECL 


where c = kb , and by defining Z=K-L Hill's formulation given by 
(1.1.5) and (1.1.7) reduces to 


Z-Z 
Ga dE ir O 


Cie ed 2) rag ort 5 


where Z = Ko L, - In the case of stimulation by a step current of 


magnitude IT, administered at t =0 , (1.5.12) is given by 


L=-.Z 
oe) 
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dz ab! 
(785913) Be arctocey & 


where 6(t) is the Dirac delta function, and has solution 
f 
} 
| 
(125 uit) d 


According to Katz" assumption, repetitive impulses occur when 
K=K-L>0O. By (1.5.14), repetitive impulses occur only if 


clan <a and only during the time interval 0 <t < ty where 


~el, 
CL. 5. 1s) t, = nae! z 


O 


is the duration of the train of impulses. This Hill-Katz model can be 
used to introduce adaptation variables into an excitation model by 
assuming that Z - Zo is the effective stimulating current instead of 


hand that <= Zo is the threshold for repetitive behavior. In applying 
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this procedure to the Hodgkin-Huxley formulation, FitzHugh (unpublished) 
succeeded in producing only a finite train of impulses to stimulation by 
step current irrespective of its magnitude. This procedure does not, 
however, resolvee the problem completely since it is impossible to 
stimulate all experimentally observed features of a finite train of 


impulses (FitzHugh: 1969). 


In 81.2 it was pointed out that the Hodgkin-Huxley equations 
are a model for the squid nerve axon with the experimentally induced 
condition of a uniform membrane potential over the entire axon membrane. 
This interpretation also applies to the models of excitation considered 
ioS1.35, S1.4 and Sie. including tne generalized model. The models of 
excitation are therefore restricted in three major ways. First, the 
formulations do not take into account the passive electrical behavior 
of the membrane to electrical stimulus. Secondly, since the squid nerve 
axon is a representative of only one class of excitable tissue, the 
models do not necessarily apply to other excitable tissue which is 
structurally different from that of the squid nerve. Finally all the 
models account for nervous stimulation in terms of an electrical 
stimulation administered to the axon. Thus the models neglect to 
account for the effect of the soma-dendritic, through which the nerve 
normal receives stimulation, on the excitable behavior of the nerve. 
The three remaining sections of Chapter I describe the attempts which 


have been made to overcome these three restrictions of the excitation 


model. 
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81.6 Mathematical Models of Action Potential Propagation 


To this point the models of nerve excitation have been 
restricted to the special case for which the membrane potential is 
assumed to be uniform over the entire axon. The models have not included 
the passive electrical properties of the nerve resulting from the 
administration of a stimulating current. Normally during the application 
of an external current at a point on the squid nerve, there are circulating 
currents which flow along and through the membrane. It is these circulating 
currents which account for the propagation over the entire nerve of an 


action potential initiated at a point on the axon. 


The distribution of the membrane potential over the membrane 
of the squid nerve axon has been modelled by assuming the axon possesses 
the properties of a core conductor. To obtain the mathematical descrip- 
tion, the axon is considered to be a cyclindrical cable-like structure of 
infinite length consisting of a conducting core, the axoplasm, which is 
separated from the extracellular conducting medium by a resistive and 


capacitative sheath, the nerve membrane. 


A number of assumptions have been placed on this basic model in 
order to simplify the mathematical formulation. First, the axoplasm and 
extracellular fluid are considered to be ohmic conductors. Secondly, by 
assuming the axon to be sufficiently thin and the axoplasmic and 


extracellular fluid resistances to be low compared with the membrane 
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resistance it is possible to ignore the radial currents in both media. 
This restriction permits the three dimensional problem to be considered 
in one dimension. Both assumptions have been shown not to introduce 
serious errors in the calculations (Hodgkin and Rushton: 1946). A 
further simplification results by considering the extracellular medium 
to be of sufficient volume in order that the extracellular fluid is 
essentially isopotential. This final assumption makes it possible to 
neglect the extracellular fluid resistance which will be negligible 


under the circumstances (Noble: 1966). 


With these assumptions, the current i, flowing through the 
axoplasm is parallel to the axon and is related to V , the displacement 
of the membrane potential from the resting potential, by Ohm's law. 


Hence 


CE26.1) —=- is 


where x is the distance along the axon, and te is the resistance of 


the axoplasm per unit length of axon. Kirchhoff's law 


(i.632) ae 


where iy is the membrane current density per unit length of axon, is 
satisfied at all points on the axon except where the external stimulus 


is applied. Equations (1.6.1) and (1.6.2) can be combined to yield 
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for which the boundary conditions are given by 
V=0O at t=0 for -~<x<o 
(1.6.4) itand 


tive t S10 Vifor) OF ec ioune 
x>too uy 


For practical application it is convenient to reformulate (1.6.3) in terms 


of more basic quantities. This is accomplished by defining 
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where d is the diameter of the axon, po is the specific resistivity of 
the axoplasm, and I is the membrane current density per unit area of 
membrane. The formulation for the distribution of the membrane potential 
over the axon is therefore given by 


d 
(1.6.6) i te aoe 


It is possible to make (1.6.6), the cable equation, more 
explicit by expressing the membrane current I in terms of its con- 
stituent parts. Since the membrane can be considered to consist of a 


capacitor and conductor in parallel, then (1.6.6) can be expanded to 


(1.6.7) didel, ay 
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where Cu is the capacity of the membrane per unit area, I, is the 
ionic current density and I. is the current density of the applied 
external current. A’ special case of (1.6.7) is that of a linear’ cable, 
for which I is a linear function Of Vv’. iii this case at it is 


assumed that the external current is applied at x = 0 , then at each 


point x #0 the membrane potential is determined by 


g 
d av ov, V 
(1.6.8) eee ea eee 
40 3x2 M ot Ry 


Where Ry is the membrane resistance per unit area of membrane. To 
solve (1.6.8) it is convenient to consider the transformation determined 


by defining 


x 
A 
(1.6.9) 4 and 
| T== 
i 
d 
where T = Row is the membrane time constant and i = aS is the 


membrane space constant. The resulting transformation of (1.6.8), given 
by 


(1.6.10) 27-2 -v=0, 


is amenable to solution by Fourier transform methods. An important 


application of cable equation (1.6.10) will be considered in §1.8. 


To model the behavior under electrical stimulation of a squid 


nerve axon, for which the membrane potential is not uniformly maintained 
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over the entire axon, it is necessary to consider not only the excitable 
properties of the nerve but also the passive electrical properties. 

Under these circumstances the mathematical formulation describing the 
axon results by combining the cable equation which determines the passive 
behavior of the nerve to electrical stimulation with a formulation des- 


cribing the excitatory behavior. 


By combining cable equation (1.6.6) with the Hodgkin-Huxley 
equations, Hodgkin and Huxley (1952d) attempted to calculate the form of 
an action potential which is propagated along the axon. This system of 
equations is equivalent to considering the Hodgkin-Huxley equations with 


(1.1.4) replaced by the partial differential equation 


2 
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(1.6.11) s 
The set of equations consisting of (1.6.11) and (1.1.5) - (1.1.12) is 
referred to as the Hodgkin-Huxley partial differential equations. To 
reduce the difficulties of solving the partial differential equation 
(1.6.11), Hodgkin and Huxley considered only the steady state propagation 
of the action potential for which the wave form and velocity are constant. 


For steady-state propagation, (1.6.6) is replaced by 


Z 

d 3Vv 

Gi 76-12) I = ae 
406 ot 


where §@ is the constant velocity of the propagated action potential 
and (1.6.11) of the Hodgkin-Huxley partial differential equations is 


replaced by the ordinary differential equation 
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Using the numerical methods of Hartree (1932-33) to solve the system of 
equations (1.6.13) and (1.1.5) - (1.1.12), Hodgkin and Huxley (1952d) 
determined that the form, amplitude and velocity of the calculated 
propagating action potential are in close agreement with experimentally 


recorded propagating action potentials. 


Further analysis of the Hodgkin-Huxley partial differential 
equation with (1.6.13) substituted for (1.6.11) has been carried out by 
Huxley (1959a,b) using a digital computer. In particular, the effect 
of temperature on the conduction velocity was calculated. Hodgkin and 
Huxley (1952d) had determined that in order for the Hodgkin-Huxley 
equations to include the effect of temperature variation on the nerve, 
each of the rate constants Ol» Bw O Bw 1 and B. must be 


multiplied by the factor 


616,14) p= 3 


where T denotes the temperature. The numerical solution of the 

system of equations (1.6.12) and (1.1.5) - (1.1.12) with the modification 
required to incorporate the effect of temperature, indicates there exists 
a critical value Ve of v such that for” > v. the propagation of an 
action potential can not occur (HUxley: 1959a). Furthermore, for yj < ve 
two solutions in addition to that a eebaschaiie to the normal propagated 


action potential were obtained (Huxley: 1959b). One of these solutions 


NK ai kes 


ae y A 


tee add evick oF (EOMGERE) ‘sea val de eborjem ‘teotvomur 9 a 
(RSdes) yateol dns aigiae | (ar, Ney) (240.1) bas (el: 10. sid: | 
bedateoins ots Yo yyriaetey ‘ity, shnsat ques ‘ierbt ads seta 
Jladeom toques dike toomsaigh weeks ok ote Inkinedty Bobics watt 


abnkaine oq yoke petsranentt ; 


Let ineregteb Talstedy slot, sau aia +6 shaylone eikcts saat? 
‘to fuo balks nood aoe (hiv T) doh best am cen id: 1) dikes 
1oetas ane)  selectdene aP oes a Segre paints: i! gekuis. td, nee) 
bus cdoighO  .fedaiboles ctw eaddiley ‘remeaehwos- ods ‘Mo 9 + . ren * 

yolxul-ntdghor sid 198 aabae ak dete, baekeeasb bat (see) be 

oy ren ait go nolishray wider tnt sete yo 13a82e@. ola otf a?) Noid: Hi 
od deity ae bis ae? “a es my ya” ey a? A asrainia rents ‘e0es ants to 

| _ Totgad ait xd bobighs, ? 


Be Was re, 


yy: 


: 4 


od4 Yo nod elon Leotvemuit silt sou dwrsqmay sed esgunwb, TF ate 4 
aoltsobtesom oti dado GED + KO Lod) bem (St.0.2). onal saupe Norineiogs 
eo eles asp esteokbus aati ac to doette: a wsameqrosat of boxkoper 
mt alana a as %, 20h, sal eae Wy, te ¥ eutey, tautaeele f 
+ (aRORE aoe no fa Katsaez0q notio— 


ve nig | ay 


mb ei ibe 


78. 


represents a low amplitude propagated action potential, while the other 
represents an infinite train of oscillations which are propagated along 

the axon. Both responses are conducted at a velocity lower than that for 
the normal propagated action potential. Neither response has been observed 


experimentally. 


FitzHugh (1969) examined the propagation of an action potential 
using the Bonhoeffer-van der Pol equation (1.5.10) as the formulation 
describing nerve excitation. The model was restricted to the case of the 


steady state propagation. Thus the required formulation is 
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The numerical solution of (1.6.15) determined by digital computer 
simulations demonstrates there exists a critical value b. of @g such 
that for 6 > %, it is impossible for an action potential to propagate. 
For each 9 < d. there are two distinct solutions of (1.6.15). One 
solution corresponds to a normal propagated action potential. The second 
solution, a low amplitude action potential propagated at a velocity lower 
than that of the normal action potential, is similar to the solution 
determined by Huxley (1959b) for the system of equations (1.6.13) and 
(1.1.5) - (1.1.12). Hence appears to play the same role in the 
Bonhoeffer-van der Pol equations as the factor W , given by (1.6.14) 


does in the Hodgkin-Huxley equations. 
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In order to be able to quantitatively describe the initiation 
of an action potential at the point of stimulation and its subsequent 
propagation along the axon, the restriction that the action potential is 
propagated at a constant velocity and waveform must be removed. There- 
fore it is necessary to obtain the solution to a formulation consisting 
of an excitation equation together with cable equation (1.6.6) or one of 
its equivalent forms. A digital computer solution to such a formulation, 
the Hodgkin-Huxley partial differential equations, has been obtained by 
Cooley and Dodge (1966). The solution was determined by approximating 


(1.6.12) with a finite difference equation. 


To obtain this approximation, a squid nerve axon of finite 
length was assumed to be represented on the one-dimensional x-axis by 
the interval [0.x 44] . This interval was then partitioned into K+l 
subintervals each of length 6x. It was supposed that over each sub- 
interval k (k = 0,1,2,3,...,K) the membrane potential is uniform and 
is determined by the Hodgkin-Huxley equations appropriately subscribed 
to distinguish individual subintervals. If the external stimulus is 
administered at the point x = 0 , then the boundary condition for 


each subinterval k (k = 1,2,3,...,K) is 


CE, 6.10) . where 


d is the axon diameter, po is the specific resistivity of the axoplasn, 
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I is the membrane current density for the mee subinterval, and Vie © 
Ma and Vey are the membrane potentials for the k - 1, fo » and 

k + 1 subintervals, respectively. To account for the external stimulus 
applied at x =0 , the boundary condition for the subinterval k=0 is 


iL Js d Onn »| 
P62 = _—_— eee | 
( By ue Ox id) (225 6x 
where fs is the stimulus current density per unit area of membrane 


and I» we Vi» d and “p are as. defined Gn (1.6.26); 


The solution to this finite difference approximation of the 
Hodgkin-Huxley partially differential equation with boundary conditions 
(1I1.6.16) and (II.6.17) is obtained by numerical integration with respect 
to t . The computed results for a propagating action potential have 
been compared with those obtained by Hodgkin and Huxley for a prapagating 
action potential with constant wave form and velocity. The two cal- 
culations agreed in regard to the peak amplitude and duration of the 
axon potential. Differences were however noted in the exponentiation 


rate in the foot of the rising phase. 


Before attempting to derive excitation and propagation models 
for excitable tissue other than the squid nerve axon,in 81.7, it is of 
mathematical interest to consider the formulation for a generalized 
propagation equation. This generalization follows from combining the 
generalized excitation formulation (1.5.9) and cable equation CL.6.6) 


The resulting formulation is 
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In order to analyze this system, Evans and Shenk (1970) have considered 


the simplified system 
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Evans and Shenk have shown that with initial conditions 


V(x,0) = $(x) 
(1.6.20) ‘ and 


(x) 


W(x,0) 


the solution of (1.6.19) is exactly the solution of the system of 


integral equation 
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Furthermore, the continuous dependence of the solutions on the initial 
values was established and solutions with initial values in a restricted 
range were shown to remain in the range for all time (Evans and Shenk: 


To70).. 
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8I.7 Mathematical Models of Excitation and Propagation 
in Myelinated Nerve 


The squid nerve axon is a representative of only one type of 
excitable tissue. For other excitable tissue, which differ structurally 
from that of the squid nerve axon, the Hodgkin-Huxley equation have been 


successfully modified to describe the excitable behavior. 


In vertebrates, the axons of certain nerves are covered except 
for small regions, called nodes of Ranvier, by a insulating myelin sheath. 
This insulating sheath restricts the excitable membrane of the axon to 
the nodes of Ranvier. Voltage-clamp analysis for a node of Ranvier of 
a myelinated Xenopus laevis nerve reveals certain similarities between 
this myelinated axon and the squid nerve axon (Frankenhaeuser: 1959, 

1960, and 1962). In both nerves the membrane current consists of a 
capacitance current and an ionic current. Furthermore, the components 
of the ionic current found in the squid nerve axon are also present in 


the myelinated toad nerve axon. 


Frankenhaeuser did however find important differences between 
the voltage-clamp data for the myelinated nerve and that for the un- 
myelinated nerve. It is these differences which account for the 
modification of the Hodgkin-Huxley equations required to make the 
formulation applicable to the myelinated nerve. In addition to the 
three components of the ionic current which are common to both nerves, 


a fourth component, a nonspecific ionic current Ty carried predominantly 


16 eavw ean ying to sy hanae 
Liowtoitss ttl h dtoliy .abeeke Mie en Wott. ins 
iad ovbd noltaepe ws but shite ‘as oeete Siv*tenr bkipe ou7 io. | 

obveutedt’ akdeap scan adidas} 63 batt bos yiiv® 


ire 


Nieaxe bexevos oTs. nevus Bete Jo Sie uly ,eeterdsivey aT 
insta oi feve galgalvenh.A we yebvass 1% tobun betes. ,enwoigat 
OF AGKe oad) to Sneteien alias kien any exyndtiaat Anmeite pet 
io. Sevres a6: bom 4. (02 abawilaws qnakyssgastoy vokvent fo vote 
nQewisad aahilweLieia aAlssase elaeved eysen perry #isqoner. 4 
Weel sisenseiaaiieyT) HORS sweet Bhtips SH? bis’ creer. be eee 

io etetetioos t0eysus onssdnam ery aboyen iyo al CORE ba ae 


eat et 


Ieenogmea add ,stomisds aya ee ohnok, fo Baw anon 


‘ 
y 


one 


| Jnseesq o#ly 976 goxe évven hlape ‘an ak bad Seiecknise staat os 36 


if 
~ 


noswied, aonceradtt eb india ‘bat tevawod tab Yeauontined inst 
=au S8d4 yo? fads sftis ores hades sit p08 _—_ gunte-ehadiad oat 
sil 20% doubts a pid sed ab 3 .owran basnableige 


Asvion Hdd 02, roMeLOD Dxb pa | 
elaiicwnessg. betes ‘at sare am intone» enon Ht 


84. 


by sodium ions, was found to be present at the node of Ranvier. Hence 
for the myelinated axon the membrane current density I per unit area 


of membrane is defined by 


fetta) UV. 
ier. 1) sips MG A torral ga era ea SE at 


for which the notation is defined in the Hodgkin-Huxley equations. The 
most significant difference between the two nerve types is that the 
instantaneous current-voltage relationship for the myelinated nerve is 
linear only in the case of the leakage current. For the sodium, potassium, 
and nonspecific currents, this relationship is most accurately described 
by the constant field equations of Goldman (1943) which in the case of 


the sodium current is 


Eur 
EF [Na], - [Na], e RT 
(1.7.2) Tue 
Na Naame Rar Ey? 
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where [Na], and [Na], are the concentrations of sodium in the 


extracellular medium and axoplasm, respectively, Pua is the sodium 


permeability, F is Faraday's constant, R is the gas constant and T 


is the absolue temperature. I. and I, are given by expressions 


similar to (I.7.2) which contain the appropriate concentrations and 


permeabilities. 


By (1.7.2), the permeabilities P Pe and Py at a given 


Na’ 


time are independent of the membrane potential. Therefore in considering 


the sodium, potassium and nonspecific ionic currents, it is more 
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appropriate to use the permeability rather than the conductance to des- 
cribe the state of the respective ionic currents. As in the case of 


the conductances, the permeabilities Py? Pa and Ey are functions 


of both time and the membrane potential. To determine the functions 


PE = P(t) and Pia = Pa St) at a given membrane potential 


Frankenhaeuser used the formulations for &, = &, (t) and fue = Bu, bt) 


of the Hodgkin-Huxley equations as a guide line. The resulting 


equations are 
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respectively, where Be and Pua are constants, and O<n,mh<1l. 


The formulation for ES = ae at a given membrane potential, which 
both satisfies the experimental data and is of the form of a Hodgkin- 


Huxley conductance equation, is 
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where ia is a constant and 0< p<1. As is the case for the 
conductances of the Hodgkin-Huxley formulation, the dependence of the 
permeabilities on the membrane potential is determined by the rate 
coefficients O, and BS fon, Lyn, lisp -..... The. formulations for ..the 
rate coefficients, as functions cf the membrane potential, were 


empirically determined and are given by 
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where A, B and C are constants and V is the displacement of the 


membrane potential from the resting potential. 


Using this formulation (1.7.1) - (1.7.6), Frankenhaeuser and 
Huxley (1964) computed the membrane currents and potential changes 
occurring during a membrane action potential. The resulting computer 
simulations were in close agreement with the experimentally recorded 


results. 
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The method of propagation of a nerve action potential in a 
myelineated nerve differs significantly from that for an unmyelinated 
nerve. For the myelinated nerve the insulating qualities of the myelin 
sheath restricts the generation of membrane current to the region of 
the nodes of Ranvier. Hence during the propagation of an action 


potential the excitation jumps from node to node. 


For the purposes of theoretical treatment, the myelinated 
axon is considered to consist of passive electrical segments, representing 
the internode lengths of myelin, separated from each other by short 
isopotential areas of excitable membrane, representing the nodes. To 
obtain the mathematical description of the propagation of an action 
potential in the myelinated nerve, Goldman and Albus (1968) represented 
a Xenopus laevis nerve by the one dimensional x-axis such that there is 
a node at x =0 and the rest are evenly spaced with internodal dis- 
tance L. By assuming the axon is symmetric with respect to the origin 
frre ne the external stimulus current is applied, only non-negative values 


of x need to be considered. 


At each node k (k = 0,1,2,3,...) the membrane potential is 
determined by equations (1.7.1) - (1.7.6) with appropriate notation to 
distinguish the individual nodes. The passive electrical segments were 
assumed to be a linear core conducting cable. In order to apply cable 
equation (1.6.11) to describe the distribution of the membrane potential 
over the internodal segments, nbareies ce0ne are required to account for 


the resistance and capacitance of the myelin sheath. The internodal 
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myelin sheath consists of a Schwann cell which has been repeatedly wrapped 
around the nerve axon. The radial resistance and capacitance per unit 


length of axon resulting from such a structure are given by 


D 
| 1 d 
Cla7.7) “ and 


respectively, where ky and ky are constants which account for 
resistance and capacitance of each layer as well as the number of layers. 
In (1.7.7), D denotes the diameter of the axon including the myelin 


while d is the diameter of the axon excluding the myelin. Hence the 


cable equation for the internode segment is 
(i778) V4 


for which«the: definition of» Vjrp'jvand «x ‘are’ given in (1.6.21). 

To complete the mathematical description of the myelinated nerve it is 
necessary to consider the boundary and continuity conditions at each 
node. For k #0 _, the boundary condition for (1.7.8) at the ee node 


is 
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where i, is the membrane current. At the origin x = 0 , the node 


corresponding to k=0 , the boundary condition for (1.7.8) is 
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Meh 20!) p= vis + 2 eget 


where i. is the externally applied stimulus current at x = 0 and 

ty is the resistance per unit length of axon of the extracellular 
medium. To ensure the continuity of the potential over the entire axon 
it was assumed that 


C7 ob.) lim V(x,t) = V,(t) LOT ee On A213 ween, 
xkL 


where V(t) is themembrane potential at the eae node. 


By assuming the ratio of the axon diameter d_ to myelinated 
axon diameter D to be constant, and the internode length L_ to be 
proportional to the myelinated axon diameter D,as is the case for the 
myelinate nerve, Goldman and Albus carried out a dimensional analysis 
of the theoretical nerve. The computed results demonstrated that the 
propagation velocity of an action potential is porportional to the 
myelinated axon diameter D , in agreement with the experimental results 


for the myelinated nerve. 


The myelinated nerve axon is not the only type of excitable 
tissue for which a modified form of Hodgkin-Huxley has been successfully 
applied. Brady and Woodbury (1957) have succeeded in modifying the 
equations to describe the excitable behavior of the ventricular fibers 
of cardiac muscle of frog. Using an entirely different formu- 


lation from that of Brady and Woodbury, Noble (1960 and 1962) was able 
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to describe the behavior of the mammalian Purkinje fibers for cardiac 


muscle. 
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81.8 Mathematical Models of the Significance of the 
Soma-Dendritic Region on the Nerves Excitability 


For the models considered in 81.1 - 81.7 nervous stimulation 
has been assumed to be an electrical stimulus administered to the axon 
of the nerve. In vivo, however, a nerve receives stimulation through 
its dendritic and soma region, from other nerves. On this basis the 
nerve can be functionally divided into the axon portion and the soma- 
dendritic portion. The models of excitation and propagation presented 
in the preceding sections are models related exclusively with the function 
associated with the axon protion of the nerve. In 81.8 two mathematical 
theories are presented which model the structural and functional 
relationship of the soma-dendritic portion of thenerve to the excitability 


of the nerve. 


In a neural system consisting of many nerves, the axon of 
one nerve will terminate on the dendritic and somatic regions of other 
nerves. The region of close contact between the axon of one nerve and 
the soma and dendrites of another is known as the synapse. The influence 
of the presynaptic nerve on the dendrites and soma of the postsynaptic 
nerve can be transmitted by either electrical or chemical means, In 
the case of a "chemical synapse", the propagation of an action potential 
along the presynaptic axon results in the release of a specific chemical 
from the presynatic membrane which produces permeability changes in the 


postsynaptic membrane at the synapse. These permeability changes 
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result in ionic currents which may alter the membrane potential in the 
neighborhood of the synapse. If the membrane potential is changed then 
this increment in potential, called the postsynaptic potential, will 
passively spread over the rest of the nerve membrane. Its distribution 
in terms of time and distance will be determined by the membrane 


resistance and capacity. (Mountcastle: 1968). 


There are two types of postsynaptic potentials each of which 
is determined by the character of the presynaptic nerve. The excitatory 
postsynaptic potential (EPSP), which is generated by an excitatory 
presynaptic nerve, is a depolarizing potential. The effect of EPSP's 
are summed both spatially and temporally. If an EPSP or the sum of a 
number of EPSP's attain an amplitude such that a critical level of 
depolarization is reached at the trigger zone of the nerve associated 
with the soma or axon, then a propagating action potential is generated. 
The other type of postsynaptic potential is an inhibitory postsynaptic 
potential (IPSP), which is produced at a synapse by an inhibitory 
presynaptic nerve. The IPSP can be either a depolarizing or hyper- 
polarizing potential depending on the value of the membrane potential. 
The effect of an IPSP is to reduce the excitability of the nerve. As 
in the case of EPSP's, the effect of IPSP's can be summed spatially and 


temporally. 


In many nerves the dendrites are branched in complex geometri- 
cal structures. Since the effect of the postsynaptic potentials decays 


with both time and distance then different spatio-temporal patterns of 
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Synaptic activity on such dendrites and the associated soma would seem 
to have different effects at the trigger zone and hence on nerve 
excitation (Rall: 1962). Two mathematical models which permit the 
prediction of the effect of various spatio-temporal patterns of synaptic 


activity have been proposed by Rall (1962, 1964). 


The experimental evidence of ionic permeability changes and 
currents at the synapse suggest that a uniform patch of dendritic 
membrane can be described by an equivalent electrical 
circuit. (Coombs, Eccles and Fatt: 1955). The model for the dendritic 
membrane proposed by Rall (1962) consists of a capacitor in parallel 
with three separate parallel conductors. The three conductors represented 
three separate ionic pathways associated with the resting membrane po- 
tential, the EPSP, and the IPSP, respectively. For this circuit the 
mathematical formulation relating the membrane potential and membrane 
current of a unit patch of dendritic membrane was proposed by Rall (1962) 
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where r, « and j represent the ionic pathway associated with the 
resting membrane, the EPSP, and the IPSP, respectively. Gos G. and > 
denote the respective separate conductances per unit area of membrane, 

and EB Ee and o denote the equilibrium potentials of the respective 
pathways. Ev is the membrane potential, Cu is the membrane capacitance 
per unit area of membrane and I is the membrane current density per 


unit area of membrane. For the model Cur Go» Gis oe Gh» E. and a 
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were assumed to be constant with respect to time and membrane potential. 


To apply (1.8.1) it was reformulated into the equivalent form given by 
(i 58.2) IRy = 3 w+ i (V - v" ) 


where 
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In order to reduce the difficulty of determining the distri- 
bution of the membrane potential over a branched dendritic tree, Rall 
(1962) considered a special class of dendritic trees which are mathe- 
matically equivalent to a cylindrical conducting core. To obtain this 
particular class, each dendritic tree is categorized in terms of the 
number n of equal branches of common radius r which are at a distance 
x from the soma. It had been assumed the cross section of the branches 
are circles, but the branches are not necessarily considered to be 


cylinders. 


To derive this classification, it is necessary to determine 
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the distribution of the membrane potential over a theoretical dendritic 
tree. The mathematics are simplified by assuming that the extracellular 
medium is isopotential and that in the intracellular medium of the 
branches, the radial current flow is sufficiently small to be neglected. 
In this case, the combined longitudinal intracellular current i, for 


d 


n dendritic branches with common radius r at a distance x from the 


soma is 
; nvr OV 
(1.8.4) i, = 4 [ 7x 
where Pq is the specific resistivity of the intracellular fluid of the 


dendrites. Furthermore, for these n branches the combined membrane 


current density I per unit area of membrane is given by 


OL 
di; dA! 
(1.8.5) I | ae “dx, 
where 
dA _ ds 
(826) aes 21rn Pie 
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and A is the membrane surface area of each of the branches. By 
differentiating (1.8.4) with respect to x and making the appropriate 
substitutions given by (1.8.5) and (1.8.6), transforms (1.8.4) to 
: ee 7 
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The change of variable given by 


oe 
E 17 Ope A 
(1.8.9) a. Sy) ‘as /* 
dx 2p, |dx! 
transforms (1.8.8) to 
1 
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If it is assumed that the conductance current of the membrane is a linear 


function of the membrane potential then 


(1.8.11) IR, = V + RyCy = 


and the combination of (1.8.10) and (1.8.11) is 
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Hence, for a given K , the distribution of membrane potential over the 
dendrites is determined by equation (1.8.12) if, by (1.8.13), its branches 


satisfy 
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where no? To and Z are the values of n, r and Z for x= 0. 


For the class of dendritic trees with K=0 , (1.8.12) reduces 
to cable equation (1.6.12) and the distribution of the membrane potential 
over the dendrites is equivalent to the distribution of the potential over 
a cylindrical core conductor whose length is expressed in units of Z. 
Furthermore, equations (1.8.6), (1.8.9) and (1.8.14) imply that in the 
case when K = 0 , equal increments of dendritic surface area correspond 
to equal increments of length of the equivalent cylinder. This fact 
permits the comparison of the result of synaptic activity in distinct 
regions of the dendritic tree which have equivalent surface area but 
are at different distances x from the soma. If each branch of the 
tree is assumed to be a cylinder, that is, as 1 in (1.8.6) then the 
class of dendritic trees with K = 0 can be extended to include den- 
dritic trees with unequal branches at a given distance x from the 
soma. In this case the constraint on the branches at a given distance 


x from the soma is 


Gin 6.5) Ny 
k 


th ; 
where Ty is the radius of the k branch and the summation over 


exhaust all branches at this value of x. 
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To determine the distribution of the membrane potential over 
a dendritic tree resulting from synaptic activity, it is necessary to 
consider both the effect of presynaptic activity on the membrane potential 
and the passive electrical properties of the membrane. Mathematically 
this is accomplished by combining (1.8.2) and (1.8.3), which account 
for the effect of presynaptic activity, with (1.8.12), which determines 
the passive distribution of the membrane potential over the dendritic 
tree. In the formulation, synaptic excitation and inhibition are 
assumed to be represented by step increases in G and C, » respectively. 
For the class of dendritic trees with K=0 _, (1.8.2) and (1.8.12) can 


be combined to yield 


2 
* 
(1.8.16) oVerw-v ) +2, 
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To investigate the structural significance of the soma- 
dendritic region Rall (1962) considered two boundary valued problems 
for the system (1.8.16) and (1.8.3). The first corresponds to the 
decay from a nonuniform membrane potential distribution over the 
dendritic tree to a uniform steady state determined by the values 
chosen for k and v" in (1.8.3). In this case the boundary 


conditions are 0 
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where L is the length of the equivalent cylinder. The solution was 
obtained by the method of separation of variables together with an 
ordinary Fourier series expansion to satisfy the initial conditions. 
The second problem examines the resulting steady-state membrane 
potential distribution over the entire dendritic tree when initially 
the tree is separated into two distinct regions of uniform membrane 
potential distribution. On the equivalent cylinder of length Z=1 , 
the separated dendritic regions correspond to the lengths of cylinder 
crom = Z' =-0 ‘to Z = Ay, and from 3:2 = As to ¥Z =L » respectively. 
Each length of the cylinder is represented by a distinct formulation of 
the system (1.8.16) and (1.8.3), defined by the regional values of k 
and v which are determined by values of Go» oe and the membrane 


constants chosen to represent the region. The boundary condition in 


this case is 


(1.8.18) ov = 0 Benz 640 anda uaa 
with both Vand ow assumed to be continuous at Z=A. The solution 


dZ 


was obtained by the method of separation of variables together with a 
generalized Fourier series expansion. In both of the problems, the 
membrane potential across the soma was assumed to be the potential 


determined for the point Z = 0 on the equivalent cylinder. 


The numerical results obtained by Rall (1962) to these two 
boundary value problems for a variety of initial conditions suggest 
there exists a functional distinction between dendritic and somatic 


symaptic activity. Dendritic synaptic activity, in particular, that 
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occurring at the dendritic periphery, seems to contribute significantly 
to the maintance and adjustment in the background level of membrane 
depolarization. Synaptic activity at or near the soma, however, seem 


to be effective in initiating precisely timed action potentials. 


The applicability of this "equivalent-cylinder model" to the 
study of the dendritic trees is seriously restricted in two respects. 
First the constraint imposed by the power of S in equation (1.8.15), 
although satisfied by the dendritic branches of motorneurons of the 
spinal cord, is not satisfied by the dendritic branches of all nerves. 
Secondly since a particular range of values of Z corresponds to a 
number of branches on the dendritic tree, it is impossible to consider 


the effect of synaptic activity which differs between branches that 


correspond to the same range of values of Z. 


To overcome these difficulties, Rall (1964) proposed a 
"compartmental model" of the soma-dendritic region. In this model the 
dendritic tree is represented by a series of discrete elements or 


compartments. 


Each isolated compartment was assumed to represent a uniform 
patch of soma-dendritic membrane. The mathematical formulation deter- 
mining the membrane potential across such a uniform patch is given by 
(1.8.1). Hence for a given compartment, denoted by Cc. » the membrane 
potential across the compartment satisfies the equivalent form of 


(1.8.1) given by 
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The notation of (1.8.19), (1.8.20) and (1.8.21) is defined as for 
equation (1.8.1) and those symbols which are not subscript with i are 
assumed to be the same for each compartment. The variable Vales 
presents the normalized value of the departure of the membrane potential 


from the resting potential such that oa 0 when Ev = E and 
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and v. = 1 when E = —£ 
al, M, € 


For the case when two compartments are directly connected 
there will be a flow of current between them. If Cc. and C, denote 


two directly connected compartments then the net flow of current from 


Cli 8.22) Bi, Vv, ~ v,) 


where Bi » the conductance from C, to C, » is assumed to be equal 


to ‘¢ » the conductance from Cc, to a - Hence, the membrane 


ae: 
potential across compartment Cc, when it is directly connected to 
another compartment, e » is determined by 

aig Bia ang 

— = (-y, - + ULV 
Saeed dt ( bat ae Wa ae V5 By” s 


aL: 5 af 


where Cc. is the capacitance of Cc. . By defining 


E.+6J, +X, 
i Ts): ot : 
1 

| 

| toh, eee 

| Cdr as 9 
(1.8.24) < Coed 

and 
eo 
ii J Cc 
L i 


then (1.8.23) is given by 
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For the general case consisting of n compartments the 


membrane potential across compartment C. LOL = JOS ee ss et 


satisfies 


CP53./26) 


where 


Choe 1) 


and 


(To6. 28) 


for two compartments C, and Ci which are not directly connected. 
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Hence for the compartment model the distribution of the membrane 


potential over the soma-dendritic region is described by the system 


of linear first order differential equations (1.8.26). 


The advantage of the compartmental model is that the region 


LOS 


of the dendritic tree which each compartment represents is not specified 


and can be arbitrarily chosen according to the need of the problem. 
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Rall (1967) solved system (1.8.26) for a chain of ten compartments of 
which the first compartment Cc, was assumed to represent the soma. 

For each of the nine compartments associated with the dendrites, various 
values of CG. and oe were chosen and the potential change produced in 


C) was determined. Transient increases of G. in any compartment 


were assumed to be of the form 


(18429) G te 


where Ty is a constant. A number of predictions regarding the effect 
of the dendrite structure on processing of presynaptic stimulation were 
determined. Three of the more important predictions are that the same 
transients applied to compartments at increasing distance from the soma 
result in potentials at the soma of decreasing size and increased rise 
time, that small conductance changes are added linearly, and that time 
course changes of the conductance transients have greater effects at 

locations on the dendrites nearer the soma. The predictions have been 


tested with only partial success in the cat spinal mororneurons (Rall 


et al: 1967). 
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CHAPTER II 


Introduction 


The nerve impulse constitutes the output signal for the nerve. 
From microelectrode recordings it is observed that in a train of 
impulses, there is variability in the intervals between successive 
intervals. The activity of the nerve therefore has a stochastic nature 
in that it is impossible to predict the exact occurrence of an impulse 
(Johannesma: 1969). Hence, it is necessary to consider the activity 
of the nerve in probabilistic terms. The experimental data must 
therefore be analyzed in terms of certain statistical quantities such 


as the distribution of interresponse intervals. 


The variability of the output signal is considered to 
import certain information regarding the generation of impulses. 
Mathematical models have therefore been developed which are capable of 
generating trains of impulse which resemble experimentally observed 
behavior. By analyzing such models it may be possible to determine the 


processes of the nerve which are responsible for such behavior. 


As pointed out above, there is a random element associated 
with impulse generation. Experimental evidence exists which suggests, 
this randomness is inherent in the synaptic input to the nerve 
(Johannesma: 1969). It is therefore of interest to determine the output 
activity of a nerve in terms of the input signal characteristics. Such 


analysis results in an input-output relationship for nerve activity which 
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can be used as the basis for the interactions between nerves of a 
mathematical nerve network model. By analyzing the behavior of small 
nerve network models for which the activity of each nerve is determined 
by this input-output relationship it may be possible to gain insight 


into the activity of more realistic nerve networks. 
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§II.1 Random Walk Model of Neuronal Variability 


As stated in the introduction to Chapter II, the variability 
in the spontaneous impulse pattern of a nerve may indicate certain 
properties of the mechanism of nerve impulse generation. In 81.1, an 
early attempt to derive a mathematical model which will generate the 


spontaneous activity of nerve is considered. 


Rodieck, Kiang and Gerstein (1962) analyzed the spike train 
recordings of the spontaneous activity of a nerve in the cochlear 
nucleus of the cat. Three measurements were considered. These are 
the interresponse time, the joint distribution of two successive 
interresponse intervals, and the distribution of time intervals between 

mt+1 ‘ : ; ‘aie : 
2 action potentials, where m is a positive integer. For each 
measurement, a histogram was plotted. These histograms give an 


estimation of the probability density function for the respective 


distributions. 


Two significant properties of the histograms enabled Gerstein 
and Mandelbrot (1964) to suggest a possible mathematical model for the 
generation of such spike train activity. First, the density function 
of the interresponse time distribution, must be able to account for 
comparatively long interresponse intervals. Furthermore, the data 
indicated for appropriate changes of the time scale, the histograms for 

; 7 : ‘ ' ; m+1 F ‘ 
the distribution of time intervals between 2 action potentials for 


each m remains approximately the same as the histogram of the inter- 
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response time distribution. For this last property, if the duration of 
successive interresponse intervals is assumed independent, then 
successive convolutions of the density function for the interresponse 
time distribution, with itself, must remain of the same functional form, 
except for linear scaling changes of the independent variable. A 
probability distribution which satisfies this criterion is said to be 


stable. 


A distribution whose density function satisfies both pro- 
perties is that of the stable distribution of order >. this sis, 
however, the probability density for the first passage time in a simple 
one-dimensional random walk which begins at a fixed distance from the 
absorbing barrier. Hence Gerstein and Mandelbrot (1964) considered 


the mathematical model which generates the spontaneous neural activity 


to be of that form. 


The random walk model considered was given the following 
physiological interpretation. The membrane potential is represented 
by a state point on the one-dimensional x-axis. The threshold 
phenomena, corresponds to the absorbing barrier for the random walk, 
and is represented by the point x =d. If the membrane potential 
reaches the threshold value, x =d _, and action potential is assumed 
to occur and the membrane potential is immediately reset to the 
resting value x= Xo: Steps taken during the random walk correspond 
to postsynaptic potentials resulting from presynaptic nerve impulses. 
An excitatory impulse produces a unit step toward the threshold,while 


an inhibitory impulse results in a unit step in the toward resting 
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potential. 


By assuming the unit step sizes to be small in comparison to 
the difference d - x and the temporal spacing between the incoming 
impulses to be small, it is possible in the limiting case to consider 
the random walk as a Wiener process (Cox and Miller: 1965). For this 
process, the stationary probability density function f(y,t,x) that 
the membrane potential is at a value y at time t after the occurrence 
of a value x , assuming the value d has not been attained in the 


time interval, satisfies the forward diffusion equation 


Z 
e151) EN ee a ea 
y ‘ 
dy 


where m_ and af are the mean and variance of the limiting process 
respectively. With the appropriate initial and boundary condition given 


by 


£(y,0,x) = S(y - x,) 


f£(-~,t,x) = 0 
CLirtee) 
and 


| £(d,t,x) = 0 


where 6 = 6(u) is the Dirac delta function; the solution (II.1.1) is 
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The probability density function for the interresponse time 
distribution is given by the probability density g(t,x,) for the 


first passage time of the absorbing barrier starting from x = x 


fo) 
Using the property of conservation of probability given by 
Sat 
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Initially, Gerstein and Mandelbrot considered the average rate 
of incoming excitatory and inhibitory impulses to be equal. For this 
case, the probability of the occurrence of an excitatory impulse and the 
probability of the occurrence of an inhibitory impulse at time t are 
equal, and m= 0 in (11.1.5). Although for this case the stability 


property holds, the temporal scaling factor required to ensure stability 
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does not closely approximate that determined from the recorded data. 

If the rate of incoming excitatory and inhibitory impulses are assumed 
tO be different, that is “m #0 in (11.1.5); then the stability 
property does not hold. However, for this case, an adequate approxi- 
mation of the variance under convolution is obtained by judicious choice 


of the parameters. 


In order to compare the theoretical predictions of the model 


Hl 
ke 


: ; ; , Z 
with the experimental data, Gerstein and Mandelbrot considered s 


and equation (II.1.5) to be of the generalized form 


Cligle6) g(t,x ) = Ke t 


where K is constant a is a constant representing the square of the 
distance between the threshold and reset value, and b is a constant 
denoting the difference in the rate of incoming excitatory and inhibitory 
impulses. For the appropriate choice of parameters it was possible to 
obtain a good fit with the experimentally determined histograms. 
Furthermore good agreement was found between the theoretical and 
experimental joint distribution of G0 successive interresponse intervals. 
Discrepancies in this case exist when the successive intervals are of 

a short duration. This was attributed to the fact that for the random 
walk model the duration of interspike intervals is independent while this 


was not the case for the experimental data. 
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Two major difficulties are inherent in the model. First it 
is impossible to obtain an exponential interresponse time distribution, 
which has been experimentally observed. Secondly, the large magnitude 
of the membrane time constant determined panne formulation has not been 
obtained experimentally for the vertebrate nervous system. Ten Hoopen 
(1966) suggests that this large value of membrane time constant can be 
attributed to the fact that each incoming impulse produces a postsynaptic 
potential which remains at a constant magnitude until it is abolished 
by an action potential. Both short comings of this random walk model 
have been overcome in the model proposed by Stein (1965), which is 


discussed.in 8II.2. 
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SII.2 Exponential Decay Model of Neuronal Variability 


The model proposed by Gerstein and Mondelbrot (1964) to 
account for the variable spontaneous behavior of nerve was based 
predominately on the mathematical properties of the experimentally 
determined distribution of interresponse times for a single nerve type. 
Although the model is successful in reproducing the distributions on 
which it is based, it is not able to generate certain other experi- 
mentally recorded interresponse time distributions. To overcome these 
difficulties it is necessary to include in such a model certain of the 
experimentally determined properties of synaptic behavior in the nerve. 
In 811.2, a model of neuronal variability proposed by Stein (1965) which 


incorporates such behavior is discussed. 


Experimental evidences indicate that postsynaptic potentials 
resulting from either excitatory or inhibitory presynaptic impulses can 
be summed, and decay approximately exponentially with time (Eccles:1957). 
Furthermore,it was found that at certain synapses, in the absence of 
neural input, excitatory transmitter is released at random and the 
resulting postsynaptic potentials decay approximately exponentially 


(Katz and Miledi: 1963). 


The model proposed by Stein (1965) incorporates these facts 
by assuming excitatory and inhibitory input impulses occur at random 


with mean frequencies n, and ns > respectively, and between impulses, 


“a rye yaa 
% ; ‘ie et " ls 
iat tare ei rn 18 7 nt 
me = : ‘- % Mf - p , if a ie t 


od (NAOLY. a6adeboge Bae nisaet9 ite baasgosy pees oh 


a 
hayad wen svaee V6 ronvetied ieroons faye aldanxee dt “a ie 
S , a 


ALAS sgn oth) th eon raga Laot Segrorloar am mm yiede 
oy svreh olgnha a 292 est) sem@gesrresar 3) coldbinitieal 
ao pooitedtqtets ofd sotovborqor at) Peteessoue at fobea ald 


" 
-Y6uxe sattia wlezdss, stessndtg) ov oite our ak af bowed ab at i 


sgan3 ainpotTeve oF ., amok: sirdt-ed etl ant SRE AN SOI babvones ' 
sii2 iv ntotwey Cshom 5 woud mt obylont os yremesaat at st a9. od 
syvten oud vt) soFvetded siaqenve to asl SeEqotg hegtoreteb elena : 
Jotdw (t80L) whats vd hawoqory Gaby dade loneweb ie — r) “a ate 


bovdvoelb, wt dobondinth i #9: 


< 


eiebiasaog sLigenvesteog téid stentha? esonmhivs ‘aelals sur 


cao saelyumt oldgenyeéxg yao tdidat 70 Utore Seu. tektites mork sake : 
- 9 


rte ae 

(Yee! ie wt) emi? datw VLhansponies ylotsaixorqys yess bre ne 
io sonsade odd mt (eomgeus wieitzeo 26 arta bnvol enw 2h ne 

r] 


ad? has wohinx te. bordatey af Isitiosoays vaopsdiose etal i 


vi labsastioyxe sical qsoub dfehstdetor shaigameereog dias ne 
] a ii 
Yh, (HORE bean baa a2) 


atopt sagtis eesatoy roan) (2b0D) ubsae. hig coh rab 
wobtay 18 Yno20 noe Laqmt suqet credtdbana | eumiahsin'% 
,sseluent souwisd bes bos ,wiewlaoagees , — ho Pi bers tie — vay oy 
a "y ve 1 
x - | 7 7 a bu Fey a 


i ical : aan MM 


" i 
ae iY win 


114. 


the depolarized subthreshold membrane potentials decays exponential 
with time constant T. Each excitatory impulse was assumed to result 
in a unit depolarization while each inhibitory impulse produces c 
units repolarization. If the threshold value of membrane potential is 
reached an action potential is assumed to occur, where upon the 


membrane potential is immediately reset to the resting level. 


With these assumptions the probability distribution F(v,t) 
that the displacement V of the membrane potential from the resting 


potential is less than or equal to the value v at time t satisfies 


oF Guster? At) oP suet) = ohh. = (n, + n, )At][F(v + Av,t) - F(,t)] 
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and 


enh Cd St Oe). Hal edt) n At [F(d,t) - F(d - 1,t)] 


where d is the difference between the threshold potential and resting 
membrane potential, and n At and n,At are the probabilities of the 
occurrence of an excitatory and an inhibitory presynaptic impulse, 
respectively, during the interval (t,t + At) . Since a depolarized 
subthreshold membrane potential decays exponentially toward the resting 


mambrane potential with time constant T then 
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The solution in (II.2.4) has not been obtained in a closed 
form. In order to investigate this formulation, two special cases of 
(11.2.4), which permit analytic treatment, have been considered. These 
cases result from assuming either that the action potential does not 
occur or that the decay of a depolarized subthreshold membrane potential 


is negligible. 


By assuming the action potential does not occur, that is 


letting d+ in (11.2.4), reduces the formulation for F(v,t) to 


(11.2.5) TRS -» See - n,[F(v,t) - F(v ~ 1,t)] 


+ n.[F( +'c,t) —- F(v,t)] for v<- 


The solution to (11.2.5) is obtained by considering the transformed 


equation 
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with boundary condition 


GEE. 2.7) C(p,t.) = 1 
where 

st iv 
(EL, 2.8) Cipeo =|) ve bapiy.e) 


is the characteristic function for the distribution F(v,t) . Equation 
(11.2.6) can be reformulated into two ordinary linear differential 


equations for which the solution with boundary condition (II.2.7) is 


given by 
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To examine the properties of this and the other special cases 
of (11.2.4) considered, Stein attempted to determine, where it is possible, 
the interresponse time distribution and mean frequency of response as 
a function of the input frequency. For the case in which the action 
potential is not assumed to occur, it is possible to examine the inter- 
response time distribution by determining the mean and variance of the 


membrane potential distribution F(v,t) . Since 
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Equations "GFh.2, 20) > Cri. 2. Fis) "and" “(Ti o2 kh) Wead "to two 
important conclusions with regard to the interresponse time distribution. 
By (11.2.10), as t + © , the mean, variance and other cumulants of the 
distribution of the membrane potential with time approach steady state 
levels. The distribution of membrane potential, therefore, becomes 
time-dependent and for either weak excitation or predominant inhibition 
in which case the steady state mean membrane potential is always less 
than the threshold value, the conditional probability that an action 
potential occurs becomes constant for times greater than several time 
constants tT. Thus, the interresponse time distribution possesses an 
exponentiar tare. by \orecelo)y and’ Chiy. 2.25). a(t) increases at a 
rate twice that of u(t) . Hence in the case of predominant inhibition 
for which H(t) <"O"" inthe bearin Gate; the conditional probability 
that an action potential occurs can increase to significant values 


before attaining its constant value. The interresponse time distribution 
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in this case will rise to a sharp peak before approaching an exponential 


aistribut ion. 


For the special case of (11.2.4) in which the decay of a 
depolarized subthreshold membrane potential is assumed to be negligible 
Stein also assumed that the effect of excitation and inhibition are 
equal to one unit of depolarization and repolarization, respectively. 
The model therefore reduces to a birth and death process or queueing 
process with an absorbing barrier at x=d. In this case system 


(II.2.4) is replaced by 


dF, (t) 
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dF (t) 


a = - nF (t) + nF, (t) 


where F(t) is the probability that membrane potential is k units, 
for k = 0,1,2,...,d , above the resting potential. The moments of 

the interresponse time distribution for this model have been determined 
by the method of probability generating functions. The mean inter- 


response time in this case is given by 
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where a is the initial displacement of the membrane potential from 


the resting potential (Bailey: 1964). Hence 


responses is 
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By (11.2.16), for predominant excitation n,n, 


the mean frequency N of 
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If, in addition to considering the decay of the subthreshold 


membrane potential to be negligible. It is also assumed that presynaptic 


inhibition is absent, then the model reduces to a Poisson process. Hence, 


the interresponse time distribution is a gamma distribution with the 
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probability density function g(t) given by 
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where to is the duration of the absolute refractory period following 
the action potential during which time incoming impulses are assumed 
to be ineffective. The mean un and variance ue of the interresponse 


time distribution for (II.2.19) are 


d 
Gri 2. 20) u, Sia s 5 Ca 
a 
and 
Z d 
Gide 2221) 0, nuit 
n 
e 


Hence the mean frequency of responses is 


(11.2. 22) eng yea 
U, tn 


For the case in which potential decay is assumed to be 
negligible and inhibition to be absent, if d is less than or equal 
to one unit of depolarization then the interresponse time distribution 
is exponential. Hence, by considering the incoming impulses to be 


randomly distributed with time, overcomes one of the difficulties 
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present in the model proposed by Gerstein and Mandelbrot (1964). 


The final special case for (11.2.4) considered by Stein 
A eenes that potential decay and the threshold are present but that 
presynaptic inhibition is not. Solutions of the formulation in this 
case were obtained by computer simulation. for each choice of the 
parameters n> d and a considered in the simulations, it was 
possible to determine parameters n> di cara to for a gamma 


distribution which would accurately approximate the simulated results. 


To determine analytically the mean frequency of responses 
for this model it is necessary to reconsider equation (1.2.12) which 
determines the mean level of the membrane potential in the absence of the 
occurrence of an action potential. In the absence of presynaptic 


inhibition (II.2.23) reduces to 
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membrane potential is sufficiently small, then the mean time U for 
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the occurrence of an action potential would correspond to the time t 


after the refractory period tS when the mean value of the membrane 


potential H(t) is equal to the threshold. Hence by (II.2.23) 
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in which case the mean frequency N of responses is given by 


(2-27) oS ea ie A 
We 


Computer simulations were found to be in agreement with (11.2.26) if 


the parameters for the simulations were chosen such that nt > Ames 


The predictions of the model, proposed by Stein, and its 
special cases have been compared, where it is possible, with existing 
experimental results. Examples of the relation between N and n, 
and n, derived, in “equations .Cii.2.17).) (LI. 2618) and GLE.2.22) ‘have 
been obtained experimentally. The prediction of a gamma interresponse 
time distribution has been found for the interresponse time histograms 
for certain nerves. Furthermore, the interresponse time histograms of 
certain nerves revealed a sharp peak followed by a long tail as had be 
predicted in the case of predominant inhibition. In general, good 
comparison was found between the predictions and the experimental 


results. 


For the model proposed by Stein (1965) the analytical predic- 
tions are restricted to the cases for which either the decay of the 
membrane potential between incoming impulses is negligible, or for 
which the action potential is assumed not to occur. Ten Hoopen (1966) 


attempted to design a model which would be more tractable than the 


ass Bee 
AVA ue ; 
ran 
oo eh esc 


Ea " ; 
i ateonteern: ~ietmeeonestyaieieianse Ke 
ae: mi 4M a 
Pie, ew be >y 


nt eas pea Agiy despsyas ab ad od Bao? forow atohislumte | 93 


rob * Ty 309 foe qeeots oye eapkgalpate edd set ergs 


act bom .1is30 ve berogata taba sas 34 anotiotherg oT 


gityeine asiw ,eldiewea. ef a: sisde. baxetenton bent, seed a2) hi 


| we 


 . bite fl, aaohisd not sgiog auld to wh aa lh Speen fete aR xo 
3 ze i i. 

Bly ae 
ved (SSS. 50) bos . (ok. S. 82d) AVL Se) anoksRepe oh navtiab, Ns i 
darugesy tent mame. 6 fb aakspl bata st ght “e st? 
ty 4) ee 


ansiposern. acho sida ial bie aid, rot tile issd'agat sotsusa enn 


Woog ae al ila supbssberg 2 wana ‘ag nt bese 
{siasmyeges ela) bos edit donbong, ond, riewaod | baer enw n08 | 


1 
§ 


eehitghiges ‘eb ay 


123. 


former, but would yield similar results in these limiting cases. The 
main difference between the two models was with regard to the decay of 
the subthreshold membrane potential between incoming impulses. For the 
model proposed by Stein, the effect of each depolarizing subthreshold 
input pulse is assumed to decay exponentially with time constant T. 
In Ten Hoopen's model, the duration of the effect of each unit of 
subthreshold potential, resulting from an input impulse, is assumed to 


be exponentially distributed with mean T. 


In order to model the nerve when both excitatory and inhibitory 
incoming impulses are present, Ten Hoopen consider two cases. If the 
subthreshold membrane potential is depolarized at time t , then an 
excitatory impulse was assumed to result in a unit depolarization, while 
an inhibitory impulse returned the membrane potential one unit towards 
the resting potential. Further, the duration of each unit of depolari- 
zation resulting from an excitatory impulse was assumed to be expo- 
nentially distributed with mean t{. For a hyperpolarized somatic 
potential at time t , an inhibitory impulse further hyperpolarized the 
potential by one unit, while an excitatory impulse returned the potential 
one unit towards the resting value. The duration of each unit of 
hyperpolarization produced by an inhibitory impulse was also assumed to 
be distributed exponentially with mean t. As in the case of the model 
proposed by Stein the arrival of excitatory and inhibitory impulses are 
each assumed to be distributed randomly with mean frequencies n, and 
n; > respectively. For this model the probability FL (t) the displace- 


ment V of the membrane potential from rest is equal to k at time t 
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after an action potential has occurred at t = 0 satisfies 
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k k+1 e i T 
+ FP. CE) [a Atlan At 
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k 
a F(t) [1 - n At ][1 - n,At - 7 Atl 


k-1 
+ FL (t) [n Atl {1 - n,At - — Atl 
rong 1.2... ,d-l “where Faq 6b) = 0, 


Grats al) At ] 
U 


OF 522 29) FC t + At) = Fea (t) [n, At] (1 - n At - 
Here Ce) pace eli nates iat 
k i e 1 


+ F,(t)[1 - n,At]}[1 - a At - CO ae] 


Koel 
T 


+B, ,(t)(1 - n,Ae][n, - At] 


for ki = -1,-2,-3,5-. 5 and 

(I1.2.30) F(t + dt) = F,(t)[1 - n,Atlin, + =JAt 
+ F (t)[1 - nat]}[1 - n,At] 
+ F(t) [n, At] [n At] 


+ F_,(t){1 - n,At}{n, + S]dt 


where Telne represents the probability of the membrane somatic 


potential returning one unit toward the resting potential during the 
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interval (t,t + At) when V(t) =k. Allowing At > 0 then equations 
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(EAO24 28) RCL 92.29) randiGhie2230)* yield 
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dF (t) l 1 
CLES 2630) eT a (n, + DF, Ct) - (n, + n, )F, (t) + (n, + DF_ Ct) 
respectively. 


The probability density function g(t) for interresponse 


time distribution is given by 
(E142. 34) g(t) = ngFy_7 (0) 


By introducing the reflecting barrier at k = -r < 0 , for which 
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close approximation of the solution of F416) and hence g(t) can 


be determined. 


A comparison with the results obtain for the model proposed 
by Stein was carried out for the limiting cases as either t+ 0 or 
d*+© , With the assumption that the action potential does not occur 


the mean value H(t) of the membrane potential at time t was obtained 


using 
(1132537) Tie peas es) 
k=-—© 

With initial condition 
Clie2¢o8) H,, (0) =O. 
Hence H(t) was determined to be 

Ae 

UE 

GEIs2. 39) H(t) = tT(n, ~ n,Q -e ) 


which corresponds to equation (I1I.2.12) for the exponential decay model 
of Stein with c= 1 and tS = 0. For the case in which decay was 
neglected and 0 < k<d _ then the two models are both equivalent a 


queueing process with absorbing barrier at k=d. 


If inhibition was neglected then the probability F(t) that 


the V(t) =k after an action potential occurs at t = 0 , satisfies 
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(1.2.40) F(t + At) = Fy, (t)[1 - nde} [(k + 12%} 
ee Oey ane ee] 
k e £ 
+ F(t) (1 - n,at}{1 - «Oy 


+ F(t) Inbtl{1 - (k - 1) 2 


POL le = Ok. 2. 64. 5aeks where 


Fact) bet) Ou: 


Allowing At > 0 , then (I1.2.40) results in the set of linear 


differential equations 
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for, eksg=)0,1,2,...,d-l where 
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F_,(t) = Fy(t) 
With the initial conditions 


F (0) = 1. 
CLi.2 42) . and 


F,.(0) a «for, Rial Zeo ee. gaou 


then an exact solution can be obtained for the set of first order 
linear differential equations (11.2.41) and for g(t) = nF 4_1(t) ; 


the density function for the interresponse time distribution. 
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In order to compare the models of Stein and Ten Hoopen when 
only excitation was present, the limiting cases were again considered. 
If the decay between impulses is negligible, the two models are 
equivalent to a Poisson process and therefore yield a gamma interres- 
ponse time distribution. For the case in which the eutan potential 
was assumed not to occur then the mean value UW, (t) of the membrane 
potential at time t for the exponentially distributed decay model of 


Ten Hoopen was determined to be 


t 


(iT. 2-43) uj(t) = nr - e Aa oes 


The same value was obtained for the exponential decay model under these 
ce: ‘ 2 ; , 
conditions. The variance 0, (t) of the membrane potential at time t 


for the exponentially distributed decay model was determined to be 


t 


(i264) O° (t) =nt(l-e ‘) 


equation (I11.2.44) can be put in the equivalent form 
4 
CEs 2545) oO (t) = 


for comparison with a. (t) given by equation (I1.2.13) for the model 


proposed by Stein in the case when ™ n, = 0 and ty = 04. Since 
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then the variance in the exponentially distributed decay model is at 
the most twice that in the case of Stein's model. This implies 
the mean interresponse interval is shorter in the case of the 


exponentially distributed decay model if the threshold is reintroduced. 


To consider the effect of this descrepancy the threshold 
and decay were reintroduced into both models. For the case in which 
only excitation was present, Stein had obtained computer simulations 
for the interresponse time distribution. In each case the distribution 
can be closely approximated by a gamma distribution. Close agreement 
between the mean and variance of these gamma distributions and the 
mean and variance determined by Ten Hoopen's model under these conditions 
was obtained if the input frequencies were corrected by a factor for 


exponentially distributed decay model. 


Thus the two models are in close agreement and the 
exponentially distributed model with its advantage of an analytical 
approach becomes a useful reference for the more realistic exponential 


decay model. 


In addition to being an approximation of the exponential 


decay model, Ten Hoopen's model has further applications. Ten Hoopen 
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suggested that the exponentially distributed decay model may be 
applicable for modelling the variable duration of post-synaptic 
potentials of incoming impulses resulting from the transmission 
properties of the synaps or the distance between the synaps and the 
soma. Furthermore, it was pointed out that nonlinear summation can 


conveniently be incorporated in the model. 


In §II.1 and 8II.2 three models have been presented which 
attempt to account for the variability in the interresponse interval. 
In 81.1 through 81.4 a variety of mathematical models for nerve 
excitation were considered. For these excitation models it was possible 
to determine a general formulation for which each of the models is a 
special case. This is also the case for the models of neuronal variability 
discussed in 8II.1 and 811.2. The generalization will be discussed in 
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8II.3 Generalized Mathematical Model for Neuronal Variability 


For the models discussed in 811.1 and 811.2, the variability 
in the interresponse intervals is attributed to the fluctuations in 
Synaptic input. Therefore, inorder to determine the mathematical 
description of the signal processing aspect of neural activities it is 
necessary to account not only for the properties of the nerve but also 
the properties of the input signal. Johannesma (1969) has presented 
a formal mathematical treatment of signal processing which is based 
on five fundamental assumptions determined both by the nerve properties 
and the input signal characteristics. The resulting formulation for 
the subthreshold behavior of the nerve to synaptic input proves to be 
the general case for which the models discussed in 81.1 and §1.2 are 


special cases. 


The first assumption is that for the nerve the subthreshold 
potential behavior is completely described by the single variable V , 
which measures the displacement of the membrane potential from the 
resting membrane potential, and for which the rate of change at time t 
is dependent only on the value of V at time t and the influence i(t) , 


of other neurons at time t . This dependence is described by 


Bates ir = £(V) + g(V) + a(t) 


where f and g are differential functions. 
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The input signal i(t) represents the membrane current 
resulting from changes in synaptic permeability. Since a nerve may 
receive synaptic input from a large number of other nerves the detailed 
structure of i(t) may be very complex. To reduce the difficulties 
of attempting to describe all the complexities of such input, i(t) is 
regarded as a stochastic process which can be characterized by its most 


prominent features, the mean 
Gesu) m(t) = (4(t)) 
and incremental variance 


4-00 
(11.3.3) s(t) = J (a(eyea(e-ty) = (aCe) Mae-t)) at 


00 


To normalize this input signal to one of zero mean and unit variance, 


i(t) is considered to be of the form 


CT, See) 1'Gt)) = me) + s(t) <4 C6) 
where 
(T1752 5) j(t) = Sty) a) : 


Hence (11.3.1) is transformed to 
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a(V,t) = £(V) + g(V) + m(t) 
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In (11.3.6) and (11.3.7) the properties of the nerve are 

determined by £(V) and g(V) and the characteristics of the input 

are determined by m(t) and s(t) . The explicit behavior of j(t) 

is unknown and it may under certain circumstances be considered as a 
stochastic carrier for which its microscopical properties are irrelevant 
to the behavior of system (11.3.6) and (11.3.7). In order to determine 
the conditions under whichthis interpretation of j(t) is valid it is 
necessary to consider three time constants which are related to the 
system (11.3.6) and (11.3.7). These are the time constant TT. of the 


statistical fluctuations of the input defined by 


co 


t ig 


(11.3.8) = Sy REE an EE GY ee 
ee ry i (i(t)ei(t-t)) - (a(t) )(a(t-1)) ; 


the time constant oF of the variation of the statistical characteristics 


of the input defined by 


(II. 3.9) tT =|< 
and the time constant e for the relaxation of the system defined by 


/ P) -l, 
(11.3.10) 2 12 G+ on] 


The second fundamental assumption of Johannesma's treatment requires 
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By this assumption, for time differences At , such that Behan Jal Fee 

the microscopical properties of j(t) are irrelevant to the behavior 
of system (11.3.6) and (11.3.7) and j(t) may be replaced by the 
white noise signal w(t) which has zero mean, unit variance, and a 
delta correlation function. Equation (11.3.6) is in this case replaced 


by the Langevin equation 


Gite 33/l2) dV = a(V,t)dt + B(V,t)dW(t) 
where 
dW(t) _ 
ae 
(w(t) ) = 0 
Chl 3 e135) 
and 


(w(t)ew(tt+t)) = 6(t) . 


Thus the replacement of j(t) by w(t) implies that for At >> oe the 
process described by (I1I.3.6) is a first order Markov process described 


by Gli. «LZ. 


In order to analyze (II.3.12) it is necessary to consider the 
transition probability density function f(y,t|x,s) the probability 
density that V has a value y at time t given that it was x at 
time s . The fact w(t) hasa delta correlation function, implies 
that the transition probability density is independent of previous history. 


Hence £(y;t)x; s) satisfies the Smoluchowski equation 
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Cit 3,14) f(y,t + At|x,s) = [f(y,t + At|z,t)£(z,t|x,s)dz 


or equivalently 


5 ny een an 
134.15.) ae fCyst x,S) = ( 1) a {A (y,t)f(y,t|x,s)} 
n=1 oy 
where 
eS io) A (y,t) = lim + fe e(y + 2, t + Atly,t)d 
oJe AY? ri ur At J y Zs y,t)dz 
t 


are the incremental moments of which the first two, the drift of V 
Mis sk 7) A, (y.t) = a(y,t) 
and dispersion of V 


(LE .3.,18) Ay(yst) = b(y,t) 
are the most significant. 


It is possible to simplify the Smoluchowski differential 
equation to a second order partial differential equation. In order to 
do this it is necessary to consider the amplitude distribution of the 
white noise w(t) . If the amplitude distribution of w(t) is assumed 
to be Gaussian then w(t) is continuous and (I1.3.12) describes a 


continuous first order Markov process. In this case 
€£133219) A (y.t) = 0 FOPLACHG>US 


in (II.3.15). Hence the third fundamental assumption of Johannesma's 
model is that w(t) is continuous. The consequence of this third 


assumption is that (11.3.15) reduces to the forward differential equation 
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(II. 3.20) = £(y,t|x,s) = - e fa(y,t)£(y,t|x,s)} 


Bea il 
fo i DCYGE) OC seks) 
oy 


where a(y,t) = A, (y,t) is the drift of V and -b(y,t) = A, (yt) is 


the dispersion of V. Furthermore £(y,t|x,s) also satisfies the 


backward diffusion equation 


Z 
(1.3.21) = £(y,t|x,s) = a(x,s) £(y,tlx,s) + 5 d(x,s)—5 £(y,tx,8) « 


2 dy 
It has been demonstrated that the drift and dispersion terms of the 
diffusion equations (I1.3.20) and (I1I.3.21) are related to the signal 


characteristics a and 8 of the fluctuation equation (I1I.3.12) 


(Stratonovich: 1963). The relationship is given by 


a(y,t) = a(y,t) a {8(y,t)}? 


Gaesaeey) / and 


b(y,t) {B(y,t)}? : 


For the stationary case, in which the mean and variance of the 
input are independent of time, the transition probability density is 
f(y,t,x) , the probability density of the transition of V from x to 
y ina time interval of duration t . Diffusion equations GELs3e 20) 


and (II.3.21) in the stationary case are replaced by 


2 
(TMv2...23) = f(y,t,x) = - S {a(y)£(y,t,x)} + oe (= b(y)£(y,t,x)} 
Bf 
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and 
e) e) yj 3° 
boca) ory, t, a) = atk). f(y,t,x) > DG)—s £(y,tsx) 5 
ot ox 2 x2 
respectively. 


In order to complete the diffusion equations for the 
stationary case it is necessary to define the initial and boundary 
conditions. For the forward diffusion equation (I1I.3.23), the obvious 


initial and boundary condition which hold are 


P{y,0.x) = (y = x) 
(25) | and 


f(-~,t,x) =O. 


To obtain a second boundary condition an assumption is required. This 
fourth fundamental assumption of the formulation establishes the 
existence of a threshold value of the membrane potential for the 
initiation of an action potential. The boundary condition for (11.2.23) 


resulting from the assumption is 
GLI. 26) f(d,t,x) = 0 


where d is the difference between the threshold potential and the 
resting membrane potential. This assumption also places a restriction 
on the transition probability density function f(y,t,x). The function 
f(y,t,x) is redefined to be the probability density that V moves 
Erom x, toy fy ini ie time interval t , assuming that V has not 


attained the value d. 


The fifth fundamental assumption is included in the formulation 
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to simplify the mathematical description of repetitive activity. This 
is accomplished by ommitting all refractory properties of the action 
potential by assuming that V is returned to a reset value X, 


immediately upon the initiation of an action potential. 


With these assumptions and formulation Johannesma was able 
to show that the formulations for the subthreshold behavior of each 
of the models discussed in 811.1 and 811.2 are special cases of the 
Smoluchowski differential equations (I1.3.15) for the stationary case 


given by 


oo 
(II. 3.27) = f(y,t,x) = J eb an A_(y)£(y.t.x) 

n=1 
where ALY) is defined by (11.3.16). The mathematical differences 
between the models is expressed in their formulations for the incremental 
moments ALY . For the five basic assumptions only the third, w(t) 
is continuous, is not satisfied by each of the models. The models which 
do not satisfy this third assumption are the exponential decay model of 
Stein (1965) and the exponentially distributed decay model of Ten Hoopen 
(1966). In both cases the input signal was assumed to be discrete in 


nature. 


For the model proposed by Gerstein and Mandelbrot (1964) 


discussed in §II.1, equation (11.3.1) is given by 


Cit 3-25) —= i(t). 
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fence by “(Il oe/), Kilss.t9), (il.3.22)% and Clivorzo) the incremental 


moments are given by 


( A, (y) = a(y) =a 
2 
A,(y) =b(y) =s 
Gie3 29) 4 
| and 
AG)» 0g for, 242.3 
| 2 


where m _ and an are defined by (11.3.2) and (11.3.3). Thus the sub- 


threshold behavior as determined by the model of Gerstein and Mandelbrot 


is described by the Smoluchowski equation 
3 3 Tunney 
CEE, 3.39) EYE Cet. y= By ime (xstry)} SAG FES taut (x, tay) 
dy 


which is equivalent to (II.1.1) of the formulation proposed by Gerstein 


and Mandelbrot. 


In order to examine the models of Stein (1965) and Ten Hoopen 
(1966) in terms of the generalized formulation it is necessary to derive 
expressions for the input when it is of a discrete nature. The 


formulation determined for the stationary discrete input i(t) is 


CLES3231) i(t) =m + sw(t) 
where 
Gilg 35332) m = ) ny cy 
k 
and 
i ) 2 
Ci, 3.33) s = ny k 
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ny is the frequency of arrival of presynaptic nerve impulses at the ee 


synaps, and cL is the modality and strength of the effect of the 


impulse at the ce synaps. 


For the exponential decay model proposed by Stein, equation 


(11.3.1) is given by 


(iT, 3,38) ee v4 ieee 


Promwlt.3./), (liedslO)s CLL. 3.5.) and @LI.3. 34) ithe incremental 


moments are 


A, (y) a(y) =m ee 


A,(y) =b(y) =s 
(II. 3.35) 


ciate” Vesyonsnal a 


and 


j 
A, (y) aa for j> 3 


where m and ge are given by (I1I.3.32) and (11.3.33). Hence the 


resulting Smoluchowski equation is 


(II. 3.36) = f(y,t,x) = ee Bd GP oe 28 Oe Yn, {£(y-c, t,x) - f(y,t,x)} 
k 


for which the integrated version is given by (11.2.5) in the formulation 


derived by Stein. 


In the case of Ten Hoopen's model for which the decay of the 
membrane potential occurs in jumps with the probability of occurrence 


exponentially distributed, only the special case where all excitatory 
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and inhibitory input impulses are of equal strength was consider. 


Hence 


WEE. Oo. o)) eas £¢C 


in (II.3.32) and (II.3.33). The incremental moments in this case are 


given by 
1 Ges = Jj PA ee Uae 
( By36) A, (y) nc + n, ¢ c) (-c) : 


cGy 2 0) and = 1.2.3...3 . where ny and n, are the frequency 
of the excitatory and imhibitory impulses, respectively. The 


Smoluchowski equation for y>0O is 


(EE.3. 39) = £ Gy $t 325) 7 = (y + c)fity + @., yx) as vity, C§x) 
T nlf (y oe mt, x) ow f(y,t,x)} 


+ n, {fly TC;t, x) = LUy,tsy 


for which the integrated forms is given by (I1.2.34) of Ten Hoopen's 
formulation. For the cases in which y= 0 and y >O analogous 


Smoluchowski equations apply. 


An extension of the model proposed by Gerstein and Mandelbrot 
(1964) was also considered by Johannesma. This model, which had first 
been proposed by Gluss (1967), resulted by including in the model of 
Gerstein and Mandelbrot (1964) the proposal that subthreshold membrane 
potentials decay exponentially with time towards the resting membrane 


potential. In this case the general equation (II.3.1) is of the form 
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Since the model also satisfies the other four assumptions of the 


peverairormulatton then’ by \CLiso. 7) scl ico Loe CLLso.22)) and (1133440) 


the incremental moments are 


A, (y) = a(y) =m - os 
A,(y) = bly) = 8° 
(ae ea 
and 
Aj) =O f0r. qo) % 


Thus the Smoluchowski equation determining the subthreshold behavior 


is 


2 
2. £(y,t,x) Sates { (m at JEG cy)! ey {s£(y,t,x)} 


dy 


Cli. 3.42) 


where m and sare given by (11.3.2) and (11.3.3). 


Having established a generalized model for neuronal variability 
based on both the properties of the nerve and the characteristics of the 
input signal, it is now of interest to determine its hehavior in terms 
of the input characteristics. In 811.4, the stationary activity of two 


special cases of the Johannesma's formulation will be discussed. 
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§I1.4 Stationary Behavior of the Perfect Integrator 
and Leaky Integrator 


In order to be able to examine a mathematical model for the 
signal processing activity of the nerve it is necessary to formulate in 
terms of the model, statistical measurements of nerve behavior which are 
amenable to experimental investigation. Such measurements include the 
distribution of interresponse time, the distribution of membrane potential, 
and the mean interresponse interval. To determine the behavior of the 
general formulation proposed by Johannesma for signal processing models, 
Johannesma (1969) analyzed the stationary activity of two special cases 
of this general formulation. These are the model proposed by Gerstein 
and Mandelbrot (1964) described by (I1.3.28), (I11I.3.29) and (I1.3.30) 
and its extension, the model proposed by Gluss (1967) discribed by 
(II1.3.40), (I11.3.41) and (I1I.3.42). The models are referred to respec- 
tively as the Perfect Integrator and Leaky Integrator, a nomenclature 
which in each case describes the membrane's effect on subthreshold inputs. 
For each model it is possible to obtain certain general conclusions 


regarding the behavior of the model in terms of the input characteristics. 


To determine the required statistical measurements, it was 
necessary to consider three probability density functions in addition 
to f(y,t,x) , defined in §11.3. These are the probability density 
e(t.x) that the first action, potential occurs at time ‘to after “Vi='x , 


the probability density h(y,t,x) that V is at a value y at time t 
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after the occurrence of the value x irrespective of whether an action 
potential has occurred in the time interval, and the probability density 
n(t,x) that an action potential occurs at time t after V =x. Two 
important special cases of these functions are g(t,x,) -and n(t,x,) 
where Xo is the reset value taken by V immediately after an action 
potential occurs. g(t,x,) is the probability density for the interval 
between action potentials, while n(t,x,) is the conditional event 
density that an action potential occurs at time t after the occurrence 


of an action potential. 


Using these probability density functions, Johannesma was able 
to derive general expressions for the stationary distribution of the 
membrane potential and the moments of the interresponse time distribution 
in terms of the coefficients a(y) and b(y) of the diffusion equation 
(I1.3.23). The stationary probability density h(y) for the membrane 


potential was determined to be 


C(y) at 
(11.4.1) h(y) = lim h(y,t) = T) (x) G— fe - xe Cl) a, 
too ai DCI 
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eCy) = ok. 
> by) 
(TT. 4.4) 


y 
C(y) = f COZ) dz; 


oid. 


aolion ne tHihiode do avidomquedie 
valanot yt klkdsdotq edd) baw Eat saul 


(seein. “Dime ( stg $38 autaill iad 26 esgno Lansogs 
Anolon Ce yedts yisaotbomnt V eo ates Pen res s9au ad ‘ay: a | 
lavcotet sit ret yehenib vst lLidtéag eet at ¢ se DS ewo50 Lf 

dneva Taeetorbuns alt ab < oie SEs <pustanes0g anni 
sone tose 803 ‘wests. J smd de apes" ‘Enksansoq moLjon ne sa @ 


slic e264 snegnesiot |. adorn vy teaes si tienda sess gad ee i 
ofa to sotiudttt2ks Ueno esas ols $03, Seip LenenGE darsaay * = Beat 
abhiudivtalh oms2 .Sasagesriesar eslt! to esata | ods ‘he Eelsneieg 9 
gokteups Baton teh add Fo Card here ide amt et en9 ada 10 er 


snsndeam af? rot. CYA eazaitel csiiadedbg| 


b ns aia 
sid 4 aie Leis j= vd m= (oe oe 
€, (ek | Shes AL) 


145. 


E(z - x) is the Heaveside unit function and T, (x) is the mean 
interresponse time. Using (11.3.22), an equivalent formulation for 
hy) can be obtained in terms of the coefficients, a(y) and B(y) 


of the fluctuation equation (I1.3.12)...In this case 


Jy) 4 aT (2) 
CLL9423) h(y). = Dx, )2 ery easy ) 5 e(z - x) ie eae 
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¥ GOD.we abas ee af ) 
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and 


y 
T(y) = f y(z)dz 


In order to determine the distribution of the interresponse times it 
is necessary to determine the probability density g(t,x,) . By the 
backward diffusion equation (11.3.24) and the property of conservation 


of probability given by 
5 d 
(11.4.5) E(t.e) = = aod) EG ytenley 


g(t,x) satisfies 
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The solution of (11.4.6), with initial and boundary conditions (II.4.7), 
exists in closed form only for the cases in which the coefficients 

a(x) and b(x) are independent of x. Therefore, to be able to 
characterize the interresponse time distribution when the coefficients 
of (11.4.6) arenot constant, Johannesma determined expressions for the 
moments TO) of the interresponse time distribution. The formulation 


obtained is 


Let Gy yey ol (2) 


Gir. 8) T (8) = nl2 J VRO)7n ee) 


dz dy He Oa 24, 
x 
where I and 8 are given in (11.4.4). Since formulations (11.4.1), 
for the stationary distribution of the membrane potential, and (11.4.8), 
for the moments of the interresponse time distribution, are derived in 
terms of the coefficients of the fluctuation equation (I1.3.12) and 
diffusion equation (I1.3.23), both expressions can be applied to the two 


models under consideration. 


For the Perfect Integrator model, the probability density 
g(t,x,) for the interval between responses can be determined explicitly 
since the coefficients a(x) and b(x) of (11.4.6) are equal to the 
constants m and re » respectively. The interresponse time distribution 


is therefore 
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(11.5429) g(t,x,) = ———— 


The mean interresponse time T, (x) for the Perfect Integator is 
Outained by solving (11.478) for a(y) =m, B(y) = s., and n= 1 . 


Hence 


d- x 


(11.4.9) | T, (x,) = eT 


and the mean frequency of response N is 


out ak ag ee 
T, (x) d axe 
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Solving (11.4.3) and (11.4.4) with o(y) =m and B(y) = s , as defined 


in the case of the Perfect Integrator, yields 


cky=sy) 
h(y) = “te {e a 
fe) 


e(y—-x2) 
as ec (yd) ae E(y = x) [1 = . }} 


(1.4.11) ¢ with 


for the stationary probability density of the membrane potential. From 
(II.4.11), Johannesma concluded that, for the Perfect Integrator, the 
stationary distribution of the membrane potential depends on the ratio 
of the mean and incremental variance of the input signal rather than 


their individual values. 


In the case of the Leaky Integrator a closed solution to 


(11.4.6), which determines the distribution of interresponse times, can 
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not be obtained since coefficient a(x) given by 
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for the Leaky Integrator model, is not independent of x. However, the 


moments of the interresponse time distribution can be determined by 
(11.4.8), with the appropriate substitutions for a and 8B. The 


resulting formulation is 
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From (I11.4.15), Johannesma concluded that if time is measured in units 


of the time constant, then the moments of the interval distribution for 


a Leaky Integrator depend only on the two combinations of system and in- 


put parameters given by X and D. For the mean interresponse time 
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T, (x) numerical computations were made. The result were formulated in 
terms of three dimensionless variables. These are N = 7 » the average 
number of responses within a time constant, M = ane the steady state 
value of the somatic potential in the absence of a threshold, and 

SS” = att By curve fitting, the relationship, in terms of these 
dimensionless variables, between the mean interresponse time T, (x) and 


the mean m and incremental variance s of the input was determined to 


be most accurately approximated by 


2 
(II. 4.16) N = ecttt8S”-aN-é 


where a, 8 , and 6 are constant. The stationary distribution h(y) 
of the membrane potential for a Leaky Integrator was also determined. 


By (11.4.3) the formulation is 
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Thus it is possible to describe stationary behavior of the 
Perfect Integrator and Leaky Integrator models in terms of the input 
signal characteristics. Of interest is the relationship determined 
between the mean interresponse time and the input signal characteristics 
for the Leaky Integrator model. This relationship, defined by (I1.4.16) 
is essentially an input-output relationship for the nerve and as such 
will form the basis of the favayaneinn equation for a nerve net model 


described in 8II.5. 
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8II.5 Mathematical Model of a Nerve Network 


In the Introduction of Chapter II it was suggested that it is 
possible to derive a mathematical model for a nerve network in terms of 
the interactions between the individual members of the network. In 
order to do this it must be possible to define the activity of each 
nerve in terms of its input signal characteristics. Therefore, a 
natural starting point for the mathematical formulation is the input- 
output relationship given by (11.4.16) which was derived for the Leaky 
Integrator Model. Since this equation would describe the behavior of 
a single member of the network, it is subscripted to distinguish the 


nerves of the net. Hence (II.4.16) is 


2 
{a ,M,+8,S,-a,N,-d, } 
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The output signal from each nerve is necessarily discrete in 
nature. Therefore it is possible to give explicit expresses for signal 


characteristics M, and si . By definition M and S. are given by 
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can be expressed as 
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where a) represents the modality and strength of the input received 
by nerve i from j , and the summation is taken over all nerves j 
in the net. Thus, the stationary activity of the network consisting 


of n nerves is given by the system 
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The basic formulation (11.5.5) has been extended to include 


the effects of the refractory period, an external stimulus, and 
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to determine the dynamic behavior of the network. To account for the 
refractory period following an action potential, Johannesma multiplied 


the righthand side of (11.5.4) by 
CEE S587) 1 - P{N, 


where Py is the duration of the refractory period in unit of the 


membrane constant T . The resulting stationary equations are 
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Cowan (1972) has extended the treatment by considering the dynamic 


behavior of the network to be described by 
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where T = =. System (11.5.9) accounts for only the interaction of 
the nerves of the network. If, however, an external input to each 
neuron is incorporated, then (11.5.9) is transformed to 
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ey contains the external input as parametric excitation. 


The description of the neural network by equation (II.5.10) 
has been derived by Cowan (1968) in a more heuristic manner. For this 
derivation the mean frequency of responses N, (t) was assumed to be 


given by the logistic function 
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the smallest current which will elicit an impulse and Ny is the 


is the current which drives the nerve at the 50% rate, i, is 


corresponding firing rate. Using the model for a patch of dendritic 
membrane derived by Rall (1962), Cowan determined the mean current 


“build up" in the cell membrane of the nerve to be 
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voltage, Cu is the membrane capacitance, B54 is the conductance 
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The heuristic assumption that the dynamic behavior of the neural network 


is given by 


| x, (t) 
(I1.5.19) sig ate ey 2 - 2 
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completes the development of the neural equation. 


Cowan (1968, 1972) also determined the properties of a specific 
case of (11.5.18) for which there is no self inhibition or excitation 
and for which the net is symmetrical. Mathematically, these criteria 


are given by 
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With these assumptions the hypothetical network resembles those present 
in the cerebral cortex and thalamus. For large external input, that 


is for large Es (II.5.19) is closely approximated by 
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which is of a form closely related to the Lotka Volterra equations for 


predator-prey-interactions. 


The properties of (11.5.21) can be studied using Gibbs 
statistical mechanics. To obtain applicable form for (11.5.20), the 


transformation 
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where q; the stationary state of the network given by 


n 


Bibs 23 = 
( ) S bark ecu ) 


is applied. Hence (11.5.21) is reduced to the equivalent form 
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is the Liapunov function for the system. Since 
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then the activity of the net is asymptotically stable and oscillations 


about the stationary states gradually damp out. 


Equation (11.5.22) can be extended to account for the 


fluctuations of the activity in the neuron by assuming 
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where dB(T) is a Langevin term with variance parameter k, . The 
Fokker-Planck equation for the probability density-—function 
E(V)sVooe++sV5 oT) associated with the Langevin equation (II.5.28) is 


given by 


CEi.5, 29) 


with equilibrium solution 
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Furthermore, since G is a sum function for each individual 
variable, then 
(11.5 31) £(v,)dv, =e dv, : 


This particular result justifies the treatment of the network as a 


Gibbs ensemble. 


A number of results were obtained by Cowan (1970) using a 


Gibbs ensemble treatment. The equilibrium density was determined to be 
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B.q, 
where B(p,q) is the Euler $-function for which p = 5 ~ and 
B,Q - 4,) 
a ; 
= ear ranean and § measures the amplitude of the fluctuation of 
activity which is given by 
2 
(x; - 4,) 
x,(1 - x,) 
CTTS 5% 33) 8 = 8, ———_—_+—_; ° 
a (x ‘4 
seed ge 
x, Cl o= oc) 


Using the equilibrium diversity, it was possible to 
determine the interval density function, amplitude density function of 
membrane potential, and the mean rate at which fluctuation in activity 
occur. With these functions, it was possible to determine certain 
general conclusions regarding the net. In particular, it was found 
that the nerves of the network spend most of their time either completely 


off or else firing at their maximum rates. 
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§II1.6 Improved Nerve Network Models 


The model proposed in 811.5 for the activity of a nerve 
network forms the basis for a number of network models which have been 
considered. Each of these models extends the original model to 
incorporate certain physiological features of the nerve which are not 
inherent in theoriginal. Mathematically these improvements are 
accomplished through modifications of the basic interaction equation 


(II.5.11). In 8II.6 three improved nerve network models are considered. 


The first modification of the model was made in order to 
account for self-inhibition in the nerves. By (I1.5.10), the stationary 
activity of the net, consisting of n members, is determined by the 
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external stimulation, and 1 given by (II.5.6), represents the 
: woh ,th ; 
influence of the j nerve on i nerve. To incorporate self- 


inhibition in the formulation, it was proposed that the dynamic behavior 


of the network is described by 
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(11.6.2) . where 


the constant B, in (11.6.1) is assumed to be equal to 1 , and Y; 
is a nonnegative constant which represents the self-inhibition of the 


t 
i . nerve cell (Oguztoreli: 1972). 


Oguztoreli (1972) has investigated the behavior of the solution 


of (11.6.2) for the case in which the initial conditions are given by 


CLI. 623) x, (0) = x, 0 where 0O < Xi og <UL ISt Orr, Lama, D Om nee 


> > 


and the network model contains infinitely many elements. Using the 
contraction-mapping principle, Oguztoreli succeeded in proving there 
exists a unique solution of (11.6.2) in an interval 0 < t <h , which 
is determined by initial condition (11.6.3). The domain of h was 
determined by the conditions for the required contraction map. It was 
also demonstrated that in the interval 0 <t<h, x, (t) is monotonic. 
For small values of Ys » weak self-inhibition, x, (t) is decreasing. 


Furthermore, the mean value of the sensitivities given by 
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sum of monotonic functions then it is of bounded variation in O0<t<h 


and therefore may exhibit random oscillations. 
The steady state equations of system (11.6.2) given by 
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with initial condition (11.6.3) have also been analyzed. The existence 
of a unique solution has been demonstrated and was determined to be 


stable if the conditions 
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hold (Leung et al: 1972). By considering small perturbation of the 
steady state solutions it was possible to analyze the linear behavior 

of the solutions of (11.6.2) for small nerve networks consisting of 

one, two, and three elements. In each case, the conditions on the 
parameters were derived to determine when the solutions of (II.6.2) 

are stable or unstable and when they exhibit damped, periodic, or growing 


oscillatory behavior (Leung et al: 1972) 
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In the case of the actual nerve, with prolonged stimulation 
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the firing rate declines and may terminate. This property of the nerve 


is known as adaptation and was previously discussed in 81.5 with regards 


to excitatory models. Stein (unpublished) suggested that the most 


appropriate manner to introduce adaptation into the interaction equation 


(11.6.2) is to add an integral term which is the convolution of an 


exponential with a factor proportional to the firing rate at a previous 


time. The resulting formulation for the activity of a network consisting 


of n nerves is 
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The behavior of the system (11.6.7) has been analyzed with 


initial conditions 


x. (0) = x, where 0O < x, <1 
= td i,o 


PieoulO). <. ‘and 
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holds. The necessary and sufficient condition in order that xy remain 
bounded within the domain of definition (0,1) for the sensitivities, 
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By the Routh-Hurwitz criterion, the real part of the roots s of 
(I1.6.20) are negative provided the coefficients bo bi» and b, and 
the determinants A, and A, are positive. If the Routh-Hurwitz 
criterion is satisfied then solutions of (I1.6.12) are stable, otherwise 
they are unstable. By the definition of the coefficients for system 
(II.6.7) and by (11.6.15), which fnsures the uniqueness of solution 

for (I1.6.12), the Routh-Hurwitz criterion is satisfied in the case 

of (I1.6.20). Hence solutions of (II.6.12) are stable. The 


discriminant of characteristic equation (I1I.6.20) is given by 
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In order to examine the behavior of (I1I.6.28) it is first 
necessary to demonstrate the existence and uniqueness of the solutions 


of the steady state equations for system (I1.6.28) given by 
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consider the Banach Space X of vector functions x(t) , which have 
continuous components x(t), x(t), x(t), x, (t) such that 
G< x, (t) < 1 for i=1,2 and t >0O. Further assume that X is 


equipped with the norm 
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result from the definition of the sensitivities. 
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T will map X into X if it can be shown that z(t) is continuous 
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z(t) is continuous for t > 0. To determine the boundary conditions 
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the Routh Hurwitz criterion can be applied to determined the stability 
of the solutions for (II.6.30) in terms of the parameters. By 


rewriting (I1.6.59) as the product of two quadratic terms given by 
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where 
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it is possible, by examining the discriminant of each term, to determine 
whether the solutions of (I1.6.59) are real or complex. Thus it is 
possible to determine the parameteric conditions for the existence and 
number of oscillatory solutions of (11.6.27). The oscillatory solutions 
are either damped or growing in nature. This results from the fact 
b, > 0 which eliminates the possibility of the existence of periodic 
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solutions. 
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of (I11.6.21), which correspond to the roots s and s Of GLI.6.4) 


are stable if and only if 


0+y+ OVA > 0 
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Similarly the solutions corresponding to the roots s. and S, are 
stable if and only if 
6@+y- OVA > 0 


(rrs6.72). 


and 
yYO+trXR+OVA>0. 


In order for oscillatory solutions to exist A > 0 in both (II.6.69) 


andes. 6./0) and. for S.| and \s 


dt 2 
(in. G7),3) muy + 6.4.2VA > VA > =v 4.6 — 27x 
and for S3 and Sy, 
(11.6. 74) -y+0+2%W>v4>-y+tO-2vr. 


Thus it is possible to determine for (I1.6.27), for the case in which 

two identical nerves are mutually inhibitory or excitatory, the 

condition on the parameters in order for the solutions to be stable 

or unstable, and to be of nonoscillatory or of damped or growing oscilla- 
tory behavior. This particular special case has been solved by numerical 


means to determine its behavior (Stein et al: 1973). 
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derived in order to include the physiological phenomena of self- 
excitation as well as self-inhibition. Both phenomena have been 
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encorporated by adding convolution integrals to the basic formulation. 


The resulting interaction equation for a nerve network containing n 
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SUMMARY 


In the thesis two basic types of mathematical models for the 
activities of the nerve have been considered. Chapter I contains those 
models which account for the recorded electrical behavior of the nerve. 
Those in Chapter II are concerned with the impulse-processing behavior 


of the nerve. 


In Chapter I a number of different models for excitation in 
the squid giant nerve axon were discussed. Each model was able to 
simulate to a degree certain aspects of the axon's behavior. It was 
demonstrated in 81.5 that it is possible to give a generalized formulation 
for excitation for which the previously discussed models are special 
case. The excitation models have also been extended. By assuming the 
axon possessed the properties of a core conducting cable the propagation 
of an action potential can be modelled. This results by combining a 
formulation for excitation with the cable equation. The Hodgkin-Huxley 
formulation for excitation in the squid nerve axon has been modified to 
model excitation and propagation in myelinated nerves which are 
structurally different from that of the squid nerve. Finally in Chapter 
I, two models were considered which demonstrated that the distance between 
the synapse and soma in electronic terms has considerable influence on 
the amplitude and time course of the membrane potentials recorded at the 


soma which result from synaptic activity. 
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In Chapter II, three models which attempt to account for 
neuronal variability have been discussed. As in the case of the 
excitation models for the squid nerve axon, it was possible to determine 
a general formulation for which these three models can be considered as 
special cases. By analysing the behavior of a special case of this 
general formulation, the Imperfect Integrator, an input-output relation- 
ship between the mean frequency of response for a nerve and the input 
signal characteristics has been determined. This relationship forms 
the basic interaction equation for a nerve network model. Three 
extensions of this basic nerve network model were considered, each one 
of which incorporates certain physiological phenomena of the nerve 
not in the original model. For one of these extensions, that which 
incorporates accommodation of the nerve in the formulation, the linear 
behavior of the solutions has been analyzed for networks containing one 


and two elements. 


The review has been restricted to the analysis of the 
activities of only the single nerve. In the case of the nerve network 
model considered the activity was determined in terms of the activity 
of the individual members. A class of models which has not been 
considered are those of a nerve network for which the net is considered 
as a continuous medium (Beurle: 1956, Griffin: 1963, and Cowan: 1972). 
For these models the active fraction of nerves as a function of time 
and space is the characteristic variable. The drawback of these models 
is that there is no connection made between the activities of the 


individual nerves and gross activity of the network. The interaction 
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equation (11.5.10) derived hopefully represents the intermediary 
between the mathematical models of activity of single nerve cells and 
those of the activity of nerve networks which are considered as a 


continuous medium. 


a 


-piavmennaad wal 
burs elias avien aigaka: to iho seated 
s a» berebianise cat nits aizeerse axon 20.3 


Pi 


re ) F ® ) it. @ eee 
*,9 op bhinestiag 


wk pale? mie 
a 
phan . a 
_ 


(G4 Tiree @ | 


ica |) ht ee Pays @y,.4 wR oD eeee 
Ne ee hale ee, Sal: 


a 


EGS 


BIBLIOGRAPHY 


AGIN, D. Some comments on the Hodgkin-Huxley equations. J. Theoret. 
Biol. ¢ 161-170. .1963. 


BALLEY, N.ieJ. The Elements of Stochastic Processes. Wiley, 1964. 


BASS, L. and McILLROY, D.K. Enzyme activities in polarized cell 
membranes. Biophys. J. 8: 99-108, 1968. 


BASS, L. and MOORE, W.J. Structural Chemistry and Molecular Biology. 
Freeman, 1968. 


BERULE, R.L. Properties of a mass of cells capable of regenerating 
pulses. Trans. Roy. Soc. (London) Ser. B240: 55-94, 1956. 


BLAIR, H.A. On the intensity-time relations for stimulation by 
electric currents. J. Gen. Physiol. 15: /09-729, 1932. 


BLAIR, H.A. On the excitation of tissue by means of condenser 
discharges. »J..Gen. Physiol. 16: 177-189, 1932. 


BLAIR, H.A. Conduction in nerve fibres. J. Gen. Physiol. 18: 
125-142 , 1934. 


BLAIR, H.A. On the relation of direct currents to condenser 
discharges as stimuli. J. Gen. Physiol. 18: 755-766. 


BRADY, A.J. and WOODBURY, J.W. Effects of sodium and potassium on 
repolarization in frog ventricular fibres. Ann. N.Y. Acad. 
Si, 109% 667-692, 1957. 


COLE, K.S., ANTOSIEWICZ, H.H. and RABINOWITZ, P. Automatic 
computation of nerve excitation. J. Soc. Ind. Appl. Math. 
= ned Re fo nD IF” bo Bs 


COLE, K.S. and MOORE, J.W. Potassium ion current in the squid giant 
axon: dynamic characteristic. Biophys. J. 1: 1-14, 1960. 


GOOMBS 5 Je Seoul COLES. .J.G. wand FATT abs The inhibitory suppression 
of reflex discharges from motoneurons. J. Physiol. 130: 
396-413, 1955. 


COOLEY, J.W. and DODGE, F.A. Jr. Digital computer solutions for 
excitation and propagation of the nerve impulse. Biophys. J. 
6: 583, 1966. 


COWAN, J.D. Stochastic Models of Neuroelectric Activity. Conference 
in Statistical Mechanics. (Ed.) S.A. Rice Univ. of Chic. Press, 
1972. 


ae | bas ih aif 


saat Py 
Jered) sb {Geek teppe ra on 


DSEL GRATE easy ot ieaitgode to ajnanele oat Bees 


ifoo bestusloy ai satiety saryaell AG (YORLELO bes od 
G2! EM Re 2b  .eydgott .eemexdmom 


Vgolols galucelom bine vitelmsid Levadeuyvsde 7.0.8 HOM bas ald : 
; ae 8a0l  nsaeext ” 


grtlerstsies 2p sidaijs> eiieo eo senza Ye ex eae at a 
DERE ReWae :O8S2 .7se (odkwod) .462 swat . ebetde (to 


wad por setuet ye 70% saokieler @nti-qilemesol.o82 #0 AH 
SECI' (QSTHROT res Loteynt fed iL. .eonetrtds” stxdoale 


ronaabags 2o-sasee KU guAgRA To doLzed toe Oat m0. | Aull ata) 
RGGt  CBL-Ter OL .Soteynt .ad.t ————- rr 


oy 


a i 
-6i .Lolaxs?- neo, .% on avrer ct es ae TA ad 
se ; 


1iawhees oS adowiegia jaedkh pliable hia hal “ 
ape ‘ 
vo soteaetoq has pekbée lo e¥esTe® "WoL, poole Ay 


28V-22y hl toseydd .o68. 1, theme, Bs Beds 
bush OY coh oe melee ree a aetna | hie , 
‘ panaideas: lik <cimibaltesias taal gee ak 309 = 
fink .bggh at a0 wigbtadtoxe evasn 4e maksasogaen ° ee je 
-ChRT SLE t€ | 


ats waite? L {3R00M bas 42.4 4 i009 
8d ge sR 2 reer Van 


oe wie ot oak cote: oe ae 


* 


136. 


COWAN, J.D. A statistical mechanics of neurons activity. In some 
mathematical questions in biology. (Ed.) M. Gerstenhaber. 
Province, R.I. Amer. Math. Soc. 1-58, 1970. 


COX, D.R. and MILLER, H.A. The Theory of Stochastic Processes. 
Methuen, 1965. 


ECCLES, J.C; The Physiology of Nerve Cells. Baltimore, Maryland. 
John Hopkins Press, 1957. 


EVANS, J. and SHENK, N. Solutions to axon equations. Biophys. J. 
10: 1090-1011,. 1970. 


FITZHUGH, R. Thresholds and plateaus in the Hodgkin-Huxley nerve 
equations. J. Gen. Physiol. 43: 867-896, 1960. 


FITZHUGH, R. Impulses and physiological states in theoretical models 
of nerve membrane. Biophys. J. 1: 455-446, 1961. 


FITZHUGH, R. A kinetic model of the conductance charges in nerve 
membrane. J. Cell Comp. Physiol. 66: (Suppl. 2) 111-117, 
19658 


FITZHUGH, R. Mathematical Models of Excitation and Propagation in 
Nerve. In: H. Schwann (Ed.) Biological Engineering. 
McGraw Hill, 1969. 


FITZHUGH, R. and ANTOSIEWICZ, H.A. Automatic compution of nerve 
excitation-detailed corrections and additions. J. Soc. Ind. 
Appl. Math. 7: 447-458, 1959. 


FRANKENHAEUSER, B. Steady state of inactivation of sodium permeability 
in myelinated nerve fibres of xenopus laevis. J. Physiol. 
148: 671-676, 1959. 


FRANKENHAEUSER, B. Quantitative description of sodium currents in 
myelinated nerve fibres of xenopus laevis. J. Physiol. 
15136:491-501, 1960. 


FRANKENHAEUSER, B. Potassium permeability in myelinated nerve fibres 
of xenopus laevis. J. Physiol. 160: 54-61, 1962. 


FRANKENHAEUSER, B. Instantaneous potassium currents in myelinated 
nerve fibres of xenopus laevis. J. Physiol. 160: 46-53, 
1962. 


FRANKENHAEUSER, B. and HODGKIN, A.L. The action of calcium on the 
electrical properties of squid axons. J. Physiol. (London) 
137: 218-244, 1957. 


FRANKENHAEUSER, B. and HUXLEY, A.F. The action potential in 
myelinated nerve fibre of xenopus laevis as computed on the 
basis of voltage clamp data. J. Physiol. (London) 171: 302-315, 


1964. 


: A 7) 
Tey De 

, ' 
; we ; mess vf 


“a 


eitee of .yelvissd anowen te wiiadhes a sidehdaes re a 


sedetiasery i CH). -saabelt rake peng ite ea i 


=, ad 
= 


_axgusoor? obvebitona® 16° (roadt ad? “~ oe en ia ig 


i 4 aaa 
baslyzeM ,sromtelal ela Sexe Go vaatet et aa 
‘ tees Read ; ; Snr 
| et Ken f 
U .egdqeks  .aeoKIRupe mora - wagers’ bas ti 
OTE hee 


avron vols uci dkehol ad? mi avendelg boa eblodaavat “a QU s " 
GAEL \oed~TAH ZED tobi’ mo! ae oe 


alsboa Laos exonits ak eezate Pert on bre Reon Kt _ HOUHS 


ovrew at asgrof> sametoubrion ofy Io Sebom otjsanks A oh gk yUHS 
EESELE AS) sigue) 200" Lobes (quod Ite -sterdmoe 


at soldngaygt? Qos @olsertset Qo etebor faok+painltsail 
ertapebiged Ladtgolaid (68) aasedo®., 4° tai ¢ 
-C3Ql LIEK mien ca Ag 

AP et 


ovis 24 Solsuqwes. of smmetisa Walt SOs gOTMA baa co e0UE, 
‘bas 2968 21 ok ses cit a onl eeoktogyae: balisieb-sol: 
| hill alcatel Be | 


(iis taunrveq muihon, 36 nae OP io oppte ehsdod «fl 
‘fotowd? Whee” ‘eurTOnEK fo ops ‘evien' Sh c 


th £594 7x09 br aed to) 
totegd? tb .akweme 


GEORGE, E.P. and JOHNSON, E.A. Solutions of the Hodgkin-Huxley 
equations for squid axon treated with tetraethylammonium and 
in potassium rich media. Australian J. Expl. Biol. Med. 
Sci. 39: 275-293, 1961. 


GERSTEIN, G.L. and MANDELBROT, B. Random walk models for the 


repetitive activity of neurons. Biophys. J. 6: 41-57, 
1964, ce 


GLUSS:),  B. A model for neuron firing with expotential decay of 
potential resulting in diffusion equations for probability 
density. Bull. Math. Biophys. 29: 233-243, 1967. 


GOLDMAN, D.E. Potential impedence and rectification in membranes. 
J. Gen. Physiol. 27: 37-60, 1943. 


GOLDMAN, D.E. A molecular structural basis for the excitation 
properties of axons. Biophys. J. 4: 167-188, 1964. 


GOLDMAN, L. and ALBUS, J.D. Computation of impulse conduction 
in myelinated fibres: theoretical basis of the velocity- 
diameter relation. Biophys. J. 8: 596-607, 1968. 


GRLPPITH, J.S. On the stability of brain-like structures. 
DLophys. Jt oe 299- 300, 100. 


HAGIWARA, S. and OOMURA, V. The critical depolarization for the 
spike in the squid axon. Japan J. Physiol. 8: 234-245, 1958. 


HARTREE, D.R. A practical method for the numerical solution of 
differential equations. Mem. Manchr. Lit. Phil. Soc. 
LRT EOT POS 238 


BILL, AgV. Extension and accommodation in nerve. Proc. Roy. Soc. 
B119: 305-355, 1936a. 


HELL, A.V. The strength duration relation for electric excitation 
of medullated nerve. Proc. Roy. Soc. B119: 440-453, 1936b. 


HODGKIN, A.L. and HUXLEY, A.F. Currents carried by sodium and 
potassium ions through the membrane of the giant axon of 
Loligo. J. Physiol. (London) 116: 449-472, 1952a. 


i es 


HODGKIN, A.L. and HUXLEY, A.F. The components of membrane conductance 
in the giant axon of Loligo. J. Physiol. (London) 116: 473-496, 


1952b. 


HODGKIN, A.L. and HUXLEY, A.F. The dual effect of membrane potential 
on sodium conductance in the giant axon of Loligo. J. Physiol. 


(London) 116: 497-506, 1952c. 


tat 


ra , pal 


sd3 Yo acotiule?. Ad BO2IHOL bam & 
baw ateoactuacom iittw baveots sone hhepea sok: 
baw .fotS .fqe3 .L metinstearA ethom dolty pref cmer ~ 


od? ot wlabem Alay mobemt 8, SOMMTIMMAM ban BRB gl 
Né-h 58 00 .eycigoke  -amexune Zo brass ae 


4 uj 

to ysoeb istisetogze datv gol fosvaa 08 Isbin A 
vailidedovg tet adolteupa wolabRRLy ok gotiluesz fatonaaoq © 
FOOL {ES-ES 2 aed ite! .[lud sepbemeb 


nsonidpent ol colsasidisods bas en fatsnot0% dy 
0CL OSE 20S tokeydd. 18a! al 


uo Joxboxe sila x03 ang lewdoustia x6luselom A, a 
eel BANA 18 20 ney anos to aatseqorg 


soitoulpos iselognh Iometgazegue? 2.0 (SUGRIA boa 1G 
-ytboavise. os Fo-ak i cekpaals fesxdt? bedank iy at 
2981 .VOaa L veydgokd = .mobkseler ss3emnkb 


V9 


eoagtours sit f-hiesd 30 yIlitdsse add nO: «Bab. 
ORI ,ANErees . we) Th veedgat 


ans 26% aorgeslusiogeh Lentstao <dT BOA.  ASUMOO han 32 
BCC! ,2AMRES 2 totaydt oh. aeqnt “nates biupe arid, a Sa 
y ei a Oar 
tb sabdowed "git dean.» thd 2%, ep too lrsnxgh., ra at 
008 sltdd ,oE .avihatM ioe sanoksaups Ist 
&-StOL “Ferree api - 


Sn 


nen 
902 oH .5o4F  .syitep BE ‘nplzetieameaas Kes ankeasiex a Va am 
OCS c22E~20€ # 
at eae Lee ie 
nolIRI tows erin epics forte peas ont rnd Wek 
REST ,CedGss 5: a 
Tid YOUMUN bnew, wt gases 
figuords. avo -auteesioq 


Parken i oan mS as 


oy — 


: eee Ly fe ete thins 


> aay! ree ace be Ly 
i ) 


188. 


HODGKIN, A.L. and HUXLEY, A.F. A quantitative description of membrane 
current and its application to conduction and excitation in 
nerve. J. Physiol. (London) 117: 500-544, 1952d. 


HODGKIN, A.L., HUXLEY, A.F. and KATZ, B. Measurement of current- 
voltage relations in the membrane of the giant axon of 
Loligo. J. Physiol. 116: 424-448, 1952. 


HODGKIN, A.L. and RUSHTON, W.A.R. The electrical constants of a 
crustacean nerve fibre. Proc. Roy. Soc. Ser. B133: 444-479, 
1956. 


HOYT, -R-C. The squid giant axon. Mathematical models. Biophys. J. 
5g599-431, 1963. 


HOVE, RC. Sodium inactivation in nerve fibers. Biophys. J. 8: 
1074-1098, 1968. 


HOOPEN, M.ten Probabilistic firing of neruons considered as a first 
passage problem. Biophys. J. 6: 435-451, 1966. 


HUXLEY, A.F. Can a nerve propagate a subthreshold disturbance? J. 
Physiol. (London) 148: 80-81, 1959a. 


HUXLEY, A.F. Ion movements during nerve activity. Ann. N.Y. Acad. 
Sci. 81: 221-246, 1959b. 


JOHANNESMA, P.I.M. Stochastic Neural Activity - A Theoretical Investi- 
gation. Ph.D. Thesis, Catholic University Nijnegen, Netherlands, 


1969" 

KATZ, B. Electric Excitation of Nerve. London: Oxford Univ. Press, 
1939" 

KATZ, B. Multiple response to constant current in frog's medullated 


nerve. J. Physiol. 88: 239-255, 1936. 


KATZ, B. and MILEDI, R. A study of spontaneous miniature potentials 
in motorneurons. J. PHysiol. 168: 389-422, 1963. 


LAPIQUE, L. J. Physiol. et pathgen 9: 620, 1907. 


LEUNG, K.V., MANGERON, D., OGUZTORELI, M.N. and STEIN, R.B. On 
a class of non-linear integro-differential equations: II. 


Analise Math., 1973. 


LEUNG, K.V., MANGERON, D., OGUZTORELI, M.N. and STEIN, R85. On 
the stability and numerical solutions of two neural models. 


Utilitas Mathematica, 1973. 


848i 


onayduemn 2o cok 
at ot eh tae aie 
BS iis cat tin ni shen? bah 


“Ita Tis to toomeidkas x 2 was. bas Aub 2 ea 


all 7 se rion ys EAL ban : 9 


le (itm teat, eid To 
S2ee. “ 


vs a 253 7h8I82ROS fies AAW ye tn 
Riimded pele ceee 908 ied m, srdty Ovxram a Be 


ie | 
MN ) 


| AA mei ; 
L ,eymeote aiebee Cok temo a Oks or C3 af 8 


f eet 


Bat _ oeprtegNs -atedt? eviem Bk) wonnertioant autboe 


bes vf 
eri? « 3b bitiahdiaeie erouxen Io gator duces os ” 
| SOBRE , LAN6SER 1 SU: a tote ria i Creates), 


Sormdtuderd biodeesivdua os sIageqor4 vase Ae 
»pePOL LOHOB 1B meena 


‘Desk Yai vata, .yalivikios svxen grlwb sjosmevom a ty a 
; deter eee 


> He ie 
1a? 


deer | eee bretnO 4 stestigenl savaedt to 6 


botetivhsa @' goed mt Goren nore oy satoques 4 qt 


eladyeser oats Links 


ms eee 


89. 


LIENARD, A. Etude des oscillations enteenues. Rev. Ge. Elec. 23: 
901-912, 1928. oa 


McILLROY, D.K. A mathematical model of the nerve impulse at the 
molecular level. Math. Biosci. 7: 313-329, 1970. 


MONNIER, A.M. L'excitation electrique destissues. Paris, 1934. 
MOUNTCASTLE, V.B. Medical Physiology. St. Louis: Mosby, 1968. 


MULLINS, L.J. An analysis of conductance changes in squid axon. J. 
Gen. Physiol. 42: 1013-1035, 1959. 


MULLINS, L.J. A single channel or a dual channel mechanism for 
nerve excitation<:: J.°Gen. Physiol. 52: 550-552, 1968. 


NACHMANSOHN, D. Chemical control of the permeability cycle in 
excitable membranes during electrical activity. Ann. N.Y. 
Acad. Sci, 137: 877-900, 1966. 


NOBLE, D. Cardiac action and pacemaker potentials based on the 
Hodgkin Huxley equations. Nature 188: 495-497, 1960. 


NOBLE, D. A modification of the Hodgkin-Huxley equations applicable 
to purkinje fibre action and pacemaker potentials. J. 
Physiol. (London) 160: 317-352, 1962. 


NOBLE, D. Applications of Hodgkin-Huxley equations to excitable 
tissue. Physiol. Rev. 46: 1-51, 1966. 


OGUZTORELI, M.N. On the neural equations of Cowan and Stein. 
Utilitas Mathematica 2: 305-317, 1972. 


RALL, W. Theory of physical properties of dendrites. Ann. N.Y. Acad. 
of Sci. 96-4: 1071-1092, 1962. 


RALL, W. Theoretical significance of dendritic trees for neural input 
output relations. In: Reiss R.F. (Ed.): Neural Theory and 
Modeling, 1964. 


RALL, W. Distinguishing theoretical synaptic potentials computed for 
different soma-dendritic distributions of synaptic input. 
J. Neurophysiol. 30: 1138-1168, 1967. 


RALL, W. BURKE, R.D., SMITH, T.G., NELSON, P.G. and FRANK, K. 
Dendritic location of synapses and possible mechanisms for 
the monosynaptic EPSP in motoneurons J. Neurophysiol. 30: 
1169-1193, 1967. 


ibSK «Set tA. 293 f thtisb (Ag; é 


Joqas eee 2 oh as 
bas exosd? Ts 


102 os 1q@en Bi, 


+. 


ses ,o9t5 a .v9d ee rm si Poo 


wit te obluq@k avger od Qo Lohan fohem Lppk amon aa ) 
OVEE (RSE-RLE 2h ‘skate atswet 

ACE ,aiakt «yeoome la ab sepintoels voisuatneata a 
Bde! »ydeor ere 2 -agatkognigtt | Teotbemt” ba: cor sre i | 


ove hiupe sf psgaeis aonny 29 atoylons oh ity ee 
ARI VREM-RIOL 48h Ledagel ie 


Tol nelnatoim \saneada ies & ro Leninads eigate ‘ bad ob mai 
2001 Stee. $2 Kode 40g, ot aaa afer ve i” 64 
ct Sioyo yJhitdsamrag eis Yo Lorine: teatmed +f af 

mh. ikvisow Leaslugoals i008 eogaccersat | 
| aoe <OOR=TTS | NEE. ed 


eX? no beeed eletritsd ug tHlemsosd bas aotios . 
‘OL NOARED. 3BOL sauzailt each eups golut 


ideotioges srehieyps yolinalnbaig bol eal’ hn om 


.L .abatttdesog! taxmpesag Bab 
Sel eSteatie 


I Luss EaRS Oy agape i 
mise bas ee ea nee tate ee ' 


a ' - 


: biroq — » wD, Tike a fp 7. ’ a. i” r ie 
pee at toy ame 
tai Ta A - ( 


Jonah oa 


nea taal 9a 


130. 


RASKEVSKY, N. Outline of a physico-mathematical theory of excitation 
and inhibition. Protoplasma 20: 42-56, 1933 


RASKEVSRKY, N. Mathematical biophysics physico-mathematical foundations 
of biology, Volume 1, 1938. 


RODIECK, R.W., KIANG, N.Y. and GERSTEIN, G.L. Some quantitative methods 
for study of sponaneous activity of single neurons. Biophys. 
We tole 59s) L962, 


RUSHTON, W.A.H. Blair's ‘Condenser theory of nerve excitation’. 
J. Gen. Physiol. 17: 481-486, 1934. 


STRATONOWICH, R.L. Topics in the theory of random noise, Vol. l 
Gordon and Breach, 1963. 


STEIN, RB. A theoretical analysis of neuronal variability. 
Biophys. J. 5: 173-194, 1965. 


STEIN, R.B. LEUNG, K.V., MANGERON, D. and OGUZTORELI, M.N. 
Improved neuronal models for studying neural networks. 
Kybernetik, 1973. 


STEIN, R.B. LEUNG, K.V., OGUZTORELI, M.N. and WILLIAMS, D.W. 
Properties of small neural networks. Kybernetik, 1973. 


TASAKI, I. and HAGIWARA, D. Demonstration of two stable potential 
states in the squid giant axon under tetraethylammonium 
chloride. J. Gen. Physiol. 40: 859-885, 1957. 


PELE dis A new interpretation of the dynamic changes of potassium 
conductance in the squid giant axon. Biophys. J. 5: 163-171, 
1965; 


VAN DER POL On relaxing oscillations. Phil. Mag. 2: 978, 1926. 
WILSON, R.H. and COWAN, J.D. Excitatory and inhibitory interactions 
in localized populations of model neurons. Biophys. J. 


12: 1-24, 1972. 


YOUNG, Y. Note on excitation theory. Psychometrika, 2: 103, 1927. 


amot tebmired miner ioad icc sonata Peta ae et ty 


ebon 3 oui sp tibbsbiieiiles axe? : re 0 eMERTARAD hers aM ‘aca Pe ae 
avdqoli .sacovaso olan Jo sail cy avoaisaogs 20 vbuse, 102, Pf 
; SORE, aii fu Masta 


\Jo2toae sured So-egondd Seenshina’? @ ‘stalh) HAW 4 
PERL TREb LED ALS icon 2180, oh 


io¥ ,galeo wohiex to excels ons * prone adel ID ITWC 


yitltceiney isoetwen to obeyiegs Leoktezoods A. ‘ae 
PORE gMOLWERT 12 .l vaudaois 


“ok OadOTSRNO Bes 0 WORERWAM’..Y. x (MULE ed 
jesus eniyvhbogse «ot. efebam levotwen. _bavorgal 


Bat la Atisiasdyt n " ve 
a> SMATLICW Ee Me _1aBROTSIDO Ved OMUTL 1848 
Cue err 9094) -wltowsod Leagan saad 44 ooktxeqoss: ei 


Isitasteq efdete . i Yo \qaliser duncan G ASAIO co Yo 
initomsLedi sna? = iene Gnaky” piopa. ond 4 a ee 1 
WROTE 2BRM +O) .Lokayd? .nod ele : 
mm Peoe uy to eae pusees ait In APE EP: won & 7 he! 
IN I6f OLE 4h SUA peer a Sherpa off” xb eon 
a cdg, 
aces ate ys, aa si ination ote ot 


‘aoesenny ctorkebetink Bias Ves BkIKs sink > 
-L .aydqotd erased Labo 2e, inet | 
tel 


is inh BPD 


reer eRe ‘S asd a aie Yair 


i =e + , a : a 
A. 


rol 1.) ee oor 
rea 4 oa whe ‘y 


Hite is 


Vy | ih 
aah " 


mi a) 
re i fy if 
yal i 
{ 7 il 
Adal Pa ans UP 
do ee ee ee 
an Ue 


